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PREFACE 


This publication is designed to support 
four implementations of the Basic FORTRAN 
IV language for the IBM System/360. The 
language described is implemented for the 
IBM System/360 Operating System, the IBM 
System/360 Disk Operating System, the IBM 
System/360 Tape Operating System, and the 
IBM System/360 Basic Programming Support 


Tape System. Differences among the lan- 
guage implementations are indicated by 
footnotes. | 


The material in this publication is 
arranged to provide a quick definition and 
syntactical reference to the Basic FORTRAN 
Iv language by means of a box format. In 


addition, sufficient text to describe each 
element and examples of possible use are 
given. 


Appendixes contain additional informa- 
tion useful in writing a FORTRAN program. 
This information consists of a table of 
source program characters, a comparison of 
the four language implementations, a list 
of FORTRAN supplied subprograms, sample 
programs, and a list of FORTRAN IV features 
and statements not available in Basic 
FORTRAN IV. 


The reader should have some knowledge of 


an existing FORTRAN language before using 
this publication. A useful source of 
information is the FORTRAN IV for 


System/360 Programmed Instruction Course, 


First Edition (August 1966) 


This publication, Form C28-6629-0, combines the specifications contained 
in the following publications and associated technical newsletters: 


IBM Operating System/360: FORTRAN IV _ (E  Level Subset), Form 


C28-6513 


IBM System/360 Basic Operating System: Specifications FORTRAN IV (16K 
Disk/Tape), Form C24-5014 and Technical Newsletters N24-5041, 


N24-5069, N21-5018 


IBM System/360 Basic Programming Support: FORTRAN IV (Tape)  Specifi- 
cations, Form C28-6504 and Technical Newsletters N28-2105 and 


N28-2168 


Although the information contained in these publications and newslet- 
ters is not obsolete, the publications themselves are no longer being 
maintained. Significant changes or additions to the specifications of 
the Basic FORTRAN IV language will be reported in subsequent revisions 
or Technical Newsletters to this publication, C28-6629-0. 


This publication was prepared for production using an IBM computer to 
update the text and to control the page and line format. Page 
impressions for photo-offset printing were obtained from an IBM 1403 
Printer using a special print chain. 


Requests for copies of IBM publications should be made to your IBM 
representative or to the IBM branch office serving your locality. 


A form is provided 
comments. 


at the back of this publication for reader's 
If the form has been removed, comments may be addressed to 
IBM Corporation, Programming Systems Publications, Department D39 
1271 Avenue of the Americas, New York, New York 10020 


© International Business Machines Corporation 1966 


Forms R29-0080 through R29-0087. This 
course is available through IBM representa- 
tives. 


Compiler restrictions and programming 
aids are contained in the programmer's 
guide for the respective system. The 
appropriate programmer's guide and this 
language publication are corequisite publi- 
cations. The programmer's guides are as 
follows: 


IBM System/360 Operating System: FORTRAN 
Iv_(E) Programmer's Guide, Form C28-6603 


IBM System/360 Disk and Tape Operating 


System: FORTR N_ IV Programmer's Guide, 
Form C24-5038 


IBM System/360 Basic Programming Support 


FORTRAN IV (Tape) Programmer's Guide, 
Form C28-6583 


References are made to information con- 
tained in the programmer's guides and in 
the following publications: 


IBM System/360 FORTRAN IV Language, Form 
C28-6515 


IBM System/360 Operating System: FORTRAN 


IV (E) Library Subprograms, Form 
C28-6596 
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INTRODUCTION 


The FORTRAN language is especially useful in writing programs for 
Scientific and engineering applications that involve mathematical compu- 
tations. Source programs written in the FORTRAN language consist of a 
set of statements constructed from the elements described in this 
publication. The FORTRAN compiler analyzes the source program state- 
ments and transforms them into an object program that is suitable for 
execution on the IBM System/360. 


The IBM System/360 Basic FORTRAN IV language is compatible with and 
encompasses the American Standards Association (ASA) Basic FORTRAN, 
including its mathematical subroutine provisions. Basic FORTRAN IV is a 
subset of FORTRAN IV, as described in the publication IBM System/360 


FORTRAN IV Language. 


The Basic FORTRAN IV language can be used with the following 
compilers: 


IBM System/360 Operating System FORTRAN IV (E) Compiler 
IBM System/360 Disk Operating System FORTRAN IV Compiler 
IBM System/360 Tape Operating System FORTRAN IV Compiler 
IBM Systen/360 Basic Programming Support Tape System FORTRAN IV 


Compiler 


All of the features and facilities in Basic FORTRAN IV also exist in 
FORTRAN IV. Equivalent results from valid programs compiled by either 
Basic FORTRAN IV or FORTRAN IV are assured by: 

1. Common data formats. 
2. Common format code routines. 
3. Common calling sequences. 
4. Common libraries. 
The following features of BaSic FORTRAN IV facilitate the writing of 


source programs and reduce the possibility of coding errors: 


1. Mixed-Mode: Expressions may consist of constants and variables, of 
the same and/or different types. 


2. Spacing Format Code: The T format code allows input/output data to 
be transferred beginning at any specified position. 


3. Literal Format Code: Apostrophes may be used to enclose literal 
data in a FORMAT Statement. 


4. The A Format Code: The A format code allows reading and writing of 
character data. 


5. Scale Factor: The scale factor allows modification of the internal 
or external representation of data. 
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10. 


11. 


12. 


Variable Attribute Control: The attributes of variables and arrays 
may now be explicitly specified in the source program. This 
facility is provided by a single explicit specification statement 
which allows a programmer to: | 


a. Explicitly type a variable as integer, real, or double preci- 
sion. 
b. Specify the dimension of arrays. 


Carriage Control: The first character of a record to be printed is 
used for carriage control. 


DOUBLE PRECISION Data Type: A third type of data is available; it 
gives greater precision than real data. 


Three Dimension Arrays: An array may have one, two or three 
dimensions. 


Six Character Variable Names: The name of a variable may contain up 
to six characters. 


Direct-Access Statements: Data records may be either read from or 
written on direct-access input/output devices in an order specified 
by the user. 


Function Subprograms may return results via the argument list. 


ELEMENTS OF THE LANGUAGE 


STATEMENTS 


Source programs consist of a set of statements from which the 
compiler generates machine instructions, constants, and storage areas. 
A given FORTRAN statement effectively performs one of three functions: 


1. Causes certain operations to be performed (e.g., add, multiply, 
branch). 

2. Specifies the nature of the data being handled. 

3. Specifies the characteristics of the source program. 


FORTRAN statements are usually composed of certain FORTRAN key words 
(see Table 1) used in conjunction with the basic elements of the 
language: constants, variables, and expressions. The five categories of 
FORTRAN statements are as follows: 


1. Arithmetic Statements: Upon execution of an arithmetic statement, 
the result of calculations performed replaces the current value of 
a designated variable or subscripted variable. 


2. Control Statements: These statements enable the user to govern the 
flow and terminate the execution of the object program. 


3. Input/Output Statements: These statements, in addition to control- 
ling input/output (1/0) devices, enable the user to transfer data 
between internal storage and an I/O medium. 


4. Specification Statements: These statements are used to declare the 
properties of variables, arrays, and subprograms (such as type and 
amount of storage reserved) and to describe the format of data on 
input or output. 


ye Subprogram Statements: These statements enable the user to name and 
define functions and subroutines. 


The basic elements of the language are discussed in this section. 
The actual FORTRAN statements in which these elements are used are 
discussed in following sections. The phrase executable statements 
refers to those statements in categories 1, 2, and 3. 


Table 1. Key Words in Basic FORTRAN IV 


Pe pe a aay Fel ge eI ae pee Cee Oe Te get ee Ey ag ee ge ee age ee he nt pee Be? ane a | 
| ABS DFLOAT FORMAT READ l 
DO FUNCTION REAL l 
| BACKSPACE DOUBLE RETURN l 
| DSIGN GO REWIND l 
| CALL GOTO l 
| COMMON END l 
| CONTINUE ENDFILE TABS SIGN | 
| EQUIVALENCE IDIM SNGL | 
EXIT IF STOP l 
| DABS EXTERNAL IFIX SUBROUTINE | 
| DBLE INTEGER | 
| DEFINE ISIGN WRITE l 
| DIM — FIND l 
| DIMENSION FLOAT PAUSE | 
Date ac en a al aa ae ae gn cle cel se ar mn cae aes sacha ve ae ae a an Sa eee tae eae Se ea ced d 
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CODING FORTRAN STATEMENTS 


The statements of a FORTRAN source program can be written ona 
Standard FORTRAN coding form, Form No. X28-7327 (see Figure 1). 
FORTRAN statements are written one to a line from columns 7 through 72. 
If a statement is too long for one line, it may be continued on aS many 
as 19 successive lines by placing any character, other than a blank or 
zero, in column 6 of each continuation line. For the first line of a 
statement, column 6 must be blank or zero. 


Columns 1 through 5 of the first line of a statement may contain a 
statement number consisting of from 1 through 5 decimal digits. Leading 
zeros in a statement number are ignored. Statement numbers may appear 
anywhere in columns 1 through 5 and may be assigned in any order; the 
value of statement numbers does not affect the order in which the 
statements are executed in a FORTRAN program. Blanks may be inserted in 
statement numbers where desired. 


Columns 73 through 80 are not significant to the FORTRAN compiler and 
may, therefore, be used for program identification, sequencing, etc. 














. X2 
IBM FORTRAN Coding Form sae 
o 
PR RAM PAGE OF 
OG PUNCHING 
NSTRUCTION CARD ELECTRO NUMBER* 
STATEMENT 5 IDENTIFICATION 
NUMBER 5 FORTRAN STATEMENT SEQUENCE 
i) { 
12 3 4 5)]6]7 8 9 10 11:12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 88 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72473 74 75 76 77 78 
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Figure 1. FORTRAN Coding Form 


Comments to explain the program may be written in columns 2 through 
80 of a line, if the letter C is placed in column 1. Comments may 
appear anywhere within the source program. They are not processed by 
the FORTRAN compiler, but are printed on the source program listing. 


The format of all FORTRAN statements requires that every word of the 
statement or expression, every name (e.g., variable name or subprogram 
name), and every constant used in a statement must be terminated by one 
of the following delimiters: 


~* . ,¥ t/fV7= () § Column 73 


In a statement of more than one line, column 73 is only considered a 
Gelimiter in the last line of the statement. As many blanks as desired 
may be written with any delimiter to improve readability. In addition, 
blanks may be inserted in key words, names, and constants; the blanks 
are ignored by the compiler. However, blanks that are inserted in 
literal data either enclosed in apostrophes or used in the H format code 
are treated as blanks within the data. The following examples show how 
blanks are treated in Basic FORTRAN IV (where b represents a blank): 


GbObbTbObb2b5 

GObTOb25 treated as GO TO 25 

GOTO25 

3b74 treated as 274 

0.b36 treated as 0.36 

*TABLDE' treated as TABLDE 
CONSTANTS 


A constant is a fixed, unvarying quantity. Three types of constants 
can be used in a FORTRAN source program: integer, real, and double 
precision. 


INTEGER CONSTANTS 


Integer Constant - a whole number written without a decimal point. 
It occupies four locations of storage. 


An integer constant may be positive, zero, or negative; if unsigned, 
it is assumed to be positive. Its magnitude must not be greater than 
the maximum and it may not contain embedded commas. 


Examples: 
Valid Integer Constants: 
0 
91 
73 
~214748 3647 
2 
Invalid Integer Constants: 
0.0 (contains a decimal point) 
27. (contains a decimal point) 
3145903612 (exceeds the allowable range) 
5, 396 (embedded comma) 
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REAL CONSTANTS 


ODED EYED EN OED GE? CE SU ED ET ee wee ae ene ED aoree a we eeewe ems ease cere ees Ee eee cere a em ees “CU ee eee EE AD OE “a TD ES EE ER A ee ee ee ee ee ee ED <oe e  S E ee ee cee  ae <Re eee aes came oTES eee soe 


Real Constant: - a number with a decimal point consisting of from 1 
through 7 significant decimal digits occupying four storage loca- 
tions. A real constant optionally may be followed by a decimal 
exponent written as the letter E followed by a signed or unsigned, 
one- or two-digit integer constant. If the decimal point is not 
present, then an E exponent specifies that the constant is real. 


Magnitude: If followed by an E decimal exponent: 0 or 16763 through 
| 1663 (i.e., approximately 1075); otherwise, it may 
consist of from 1 through 7 decimal digits. 


FE ED ASE SE ANE EE ET EES OIE ARNE END SERED GUND -coNS SoTES RE cenee eoESO SEE SEER ERED CEen EERE EE ARID SRE ER ee QED CEE Se ETE SEES? EEE -HESES HN: SEED ED CTE AEE SERRE SERED GR RED SEED AED ED ER eee GEE SEE eee ome QS ARES OED sO eee ce ee acme aE GREE SED ene coed a cee ques Soe 


A real constant may be positive, zero, or negative (if unsigned, it 
is assumed to be positive) and must be of the allowable magnitude. It 
may not contain embedded commas. The decimal exponent permits the 
expression of a real constant as the product of a real constant times 10 
raised to a desired power. 


Examples: 


Valid Real Constants: 


+0. 

-999. 9999 

0.0 

5764.1 

7.0E+0 (i.e., 7.0 x 10° = 7.0) | 
19761. 25E+1 (1.e., 19761.25 x 101 = 197612.5) 
TE3 

7.0E3 (i.e., 7.0 x 103 = 7000.0) 
7.0E03 

7. 0E+03 

7.0E-03 (i.e., 7.0 x 1073 = 0.007) 


Invalid Real Constants: 


0 (missing a decimal point) 
3,471.1 (embedded comma) 
1.E (missing a one- or two-digit integer 


constant following the E. Note that it is not 
interpreted as 1.0 x 10°) 


1.2E+113 _ (E is followed by a 3 digit 
| integer constant) 
23.5E+97 (value exceeds the magnitude permitted; that is, 
23.5 x 109751663) 
~-91437.143 (exceeds the number of significant decimal 
digits permitted) 
1.0000000 (exceeds the number of significant decimal 


digits permitted) 
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DOUBLE PRECISION CONSTANTS 


|Double-Precision Constant - a number with a decimal point optionally 
{followed by a decimal exponent. This exponent may be written as the 
jJletter D followed by a signed or unsigned, one- or two-digit integer 
|constant. If the decimal point is not present, then a D_ exponent 
[specifies that the constant is double precision. A double-precision 
| constant may assume one of two forms (both forms occupy eight storage 
| locations) : 


1. From i through 7 decimal digits followed by a D_ decimal 


2. From 8 through 16 significant decimal digits optionally followed 
by a D decimal exponent. 


Magnitude: (either form) 0 or 1676? through 16°? (i1.e., approximate- 
ly 1075). 


Fo ane oo oe ee ee ee ee ee eee 


A double-precision constant may be positive, zero, or negative (if 
unsigned, it is assumed to be positive) and must be of the allowable 
magnitude. It may not contain embedded commas. The decimal exponent D, 
Similar to that for real constants (i.e., E), permits the expression of 
a double-precision constant as the product of a double-precision 
constant times 10 raised to a desired power. Note that a double- 
precision constant has more than twice the number of significant digits 
as that for a real constant, but its maximum magnitude remains the same. 


Examples: 


Valid Double-Precision Constants: 


21.98753829457168 

1.0000000 

79D3 

7.9D03 

7.9D+3 (i.e., 7.9x103 = 7900.0) 

7.9D+03 

7.9D-03 (i.e., 7.9x1073 = 0.0079) 

7.9D0 (i.e., 7.9x10° = 7.9) 

0.0D0 (i.e., 0.0x10° = 0.0) 

Invalid Double-Precision Constants 

7.9E3 (should have decimal exponent D, not E) 

7.9D (missing a one- or two-digit integer 
constant following the D) 

7.987143 (decimal exponent D must follow when less 
than 8 significant decimal digits are used) 

21.3D90 (exceeds given magnitude, i.e., 


21. 3x1099>51663) 
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VARIABLES 


A FORTRAN variable is a symbolic representation of a quantity that is 
aSSigned a value. The value may either be unchanged (i.e., constant) or 
may change either for different executions of a program or at different 
stages within the program. For example, in the statement: 


A = 5.0+B 


both A and B are variables. The value of B is determined by some 
previous statement and may change from time to time. The value of A 
varies whenever this computation is performed with a new value for B. 


VARIABLE NAMES 


ie a eet eg oe as eS Peres ee le eg eg thee eg gee ad ihe ome ah ge ee eee ete Ee Grey keyg oe 


| Definition | 


| Variable Name - from 1 through 6 alphameric (i.e., numeric, 0 | 
| through 9, or alphabetic, A through Z and $ ) characters, the first | 
| of which must be alphabetic. | 
L 


Variable names are symbols used to distinguish one variable from 
another. A name may be used in a source program in one and only one way 
(e.g., the name of a variable and that of a subprogram may not be 
identical in the same source program). A variable name may not contain 
special characters other than the blank (see Appendix A). 


The use of meaningful variable names can serve as an aid in 
documenting a program. That is, someone other than the programmer may 
look at the program and understand its function. For example, to 
compute the distance a car traveled in a certain amount of time at a 
given rate of speed, the following statement could have been written: 


X= Y¥ * Z 


where * designates multiplication. However, it would be more meaningful 
to someone reading this statement if the programmer had written: 


DIST = RATE * TIME 


Examples: 


Valid Variable Names: 
B292 
RATE 
SQ704 
$VAR 


Invalid Variable Names: 


B292704 (contains more than six characters) 
4WARRAY (first character is not alphabetic) 
SI.X (contains a special character) 
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VARIABLE TYPES AND LENGTHS 


The type of a variable corresponds to the type of data the variable 
represents. Thus, an integer variable represents integer data, a real 
variable represents real data, etc. 


The number of storage locations reserved for variables depends on the 
type of the variable. Integer ‘and real variables have four storage 
locations reserved; double precision variables have eight storage 
locations reserved. 


The ways a programmer may declare the type of a variable are by use 
of the: 


1. Predefined specification contained in the FORTRAN language. 


2. Explicit specification statements. 


TYPE DECLARATION BY THE PREDEFINED SPECIFICATION 


The predefined specification is a convention used to specify 
variables as integer or real as follows: 


1. If the first character of the variable name is I, J, K, L, M, or N, 
the variable is integer. 


2. If the first character of the variable name is any other alphabetic 
character, the variable is real. 


This convention is the traditional FORTRAN method of implicitly 
specifying the type of a variable as being either integer or real. In 
all examples that follow in this publication, it is presumed that this 
specification holds unless otherwise noted. 


TYPE DECLARATION BY EXPLICIT SPECIFICATION STATEMENTS 


Explicit specification statements differ from the first way of 
specifying the type of a variable, in that an explicit specification 
Statement declares the type of a particular variable by its name rather 
than as a group of variables beginning with a particular character. 


For example, assume that an Explicit specification statement declared 
that the variable named ITEM is real. Then ITEM is treated aS a real 
variable but all other variables beginning with the character I are 
treated as integer variables. 


These statements are discussed in greater detail in the section, 
"Specification Statements." 
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EXPRESSIONS 


Expressions in their simplest form consist of a single constant or 
variable. They may also designate a computation between two or more 
constants and/or variables. Expressions may appear in arithmetic 
statements and in certain control statements. 


Basic FORTRAN IV provides only one kind of expression: the arithmetic 
expression. The value of an arithmetic expression is always a number 
whose type is integer, real, or double precision. 


ARITHMETIC EXPRESSIONS 


The simplest arithmetic expression consists of a single constant, 
variable, or subscripted variable (see the discussion of arrays). The 
constant or variable may be one of the following types: 


1. Integer 
2.2 Real 
3. Double Precision 


If the constant, variable, or subscripted variable is of the type 
integer, the expression is in the integer mode. If it is of the type 
real, the expression is in the real mode, etc. 


Examples: 

Expression Type of Quantity Mode of Expression 
3 Integer Constant Integer 

I Integer Variable Integer 

3.0 Real Constant Real 

A Real Variable Real 

3.14D3 Double-precision Constant Double precision 
B(2*T) Double-precision Subscripted Double precision 


Variable (Specified as 
Such in an Explicit 
Specification statement) 


In the expression B(2*#I), the subscript (2*1I), which must always 
represent an integer, does not affect the mode of the expression. That 
is, the mode of the expression is determined solely by the type of 
constant, variable, or subscripted variable appearing in that expres- 
sion. 


More complicated arithmetic expressions containing two or more 
constants and/or variables may be formed by using arithmetic operators 
that express the computation(s) to be performed. 


Arithmetic Operators 


The arithmetic operators are as follows: 


Arithmetic Operator Definition 

* + Exponentiation 
* Multiplication 
/ Division 

+ Addition 

_ Subtraction 


RULES FOR CONSTRUCTING ARITHMETIC EXPRESSIONS: The following are the 
rules for constructing arithmetic expressions that contain arithmetic 
operators: 


a es 


All desired computations must be specified explicitly. That is, if 
more than one constant, variable, subscripted variable, or function 
reference (see the section "SUBPROGRAMS") appears in an arithmetic 
expression, they must be separated from one another by an arithme- 
tic operator. For example, the two variables A and B will not be 
multiplied if written: 


AxB or AB or AeB or A(B) 


If multiplication is desired, then the expression must be written 
as follows: 


A*B or B*¥A 


No two arithmetic operators may appear in sequence in the same 
expression. For example, the fo'’iowing expressions are invalid: 


A*/B and A*-B 
The expression A*-B could be written correctly as follows: 
A* (-B) 


In effect, -B will be evaluated first and then A will be multiplied 
by the result. (For further uses of parentheses, see Rule 6.) | 


The mode of an arithmetic expression is determined by the type of 
the operands (where an operand is a variable, constant, function 
reference, or another expression) in the expression. Table 2 
indicates how the mode of an expression that contains operands of 
different types may be determined using the operators: +, -, *, /. 


Table 2. Determining the Mode of an Expression Containing Operands 
of Different Types 





aaa ala a ae cana a fo ee 1 
l | | DOUBLE | 
|j+ - * / INTEGER | REAL | PRECISION | 
[INTEGER | Integer | Real | Double | 
| | | | Precision | 
~--------}------~----}-----------}-----------4 
| REAL | Real | Real | Double | 
| | | Precision | 
le ee a a a a ae NI as ———————-—-——— 
| DOUBLE | Double | Double | Double | 
| PRECISION} Precision | Precision | Precision | 
[ RRM CS mee etic ea ET etme em einer £5 can sa en a Serle ee ee 4 


From Table 2 it can be seen that there is a hierarchy of type that 
determines the mode of an expression. For example, double- 
precision data when combined with any other types of constants and 
variables results in double-precision. 


Assume that the type of the following variables has been specified 
as follows: 


Variable Names Type 


ROOT, E Real variables 
A, I, F Integer variables 
CD Double-precision variables 


Elements of the Language 15 


16 


Then the following examples illustrate how constants and variables 
of different types may be combined using the arithmetic operators: 


+, ~_ Fe *F: 


Expression 
ROOT*5 


At3 
C+#2.9D10 
E/F+t#19 
C-18.7E05 
A/I-D 
C/D 


The arithmetic operator denoting exponentiation 


Mode of Expression 
Real 


Integer 

Double Precision 
Real 

Double Precision 
Double Precision 
Double Precision 


(i.e.,**) may be 


used to combine any types of operands as shown in Table 3. 


Table 3. Valid Combinations 


Operator, 





| 

| 

| 

| Integer 
| 

| Double Precision 
L 


* & 


Exponent 


| 

| 

| 

Integer | 

r+ | Real | 
| 


Double Precision 


D 
RRND AD ERED <ENNED ARSED GEG GRE GENTE CANS EGSS -<GERED SEUED ORRED QED SEE TRS QUT NEN GUNES ~EEUTD SUEDE PCIE QOD TRAE GESED <OUTD AUTRES CONRAD CEREED END EES ER GERED SEED IED LED ED ES A eR ce 


with Respect to 


the Arithmetic 


Assume that the type of the following variables has been specified 


as follows: 


Variable Names 
ROOT,E 

Ac a. F 

Cc, D 


Type 
Real variable 
Integer variables 


Double-Precision variables 


Then the following examples illustrate how constants and variables 
of different types may be combined using the arithmetic operator, 


** , 


Examples: 


Expression 
ROOT#* (A+2) 
ROOT * * I 
[**F 
7.98E21**ROOT 
ROOT*#*2.,1E5 
A*¥#*¥E 

C¥*A 

E**C 

L*¥*C 

D**E 

D*¥*C 


Type 


(Real**Integer) 

(Real** Integer) 

(Integer**Integer) 

(Real**Real) 

(Real**Real) 

(Integer**Real) 

(Double Precision** Integer) 

(Real**Double Precision) 

(Integer**Double Precision) 

(Double Precision**Real ) 

(Double Precision**Double 
Precision) 


Result 


(Real) 

(Real) 

(Integer) 

(Real) 

(Real) 

(Real) 

(Double Precision) 
(Double Precision) 
(Double Precision) 
(Double Precision) 
(Double Precision) 


Order of Computation: Where parentheses are omitted, or where the 
expression is enclosed within a single pair of 


entire arithmetic 


parentheses, effectively the order 


performed is as follows: 


in which the operations are 


Operation Hierarchy 


Evaluation of Functions (see the ist (highest) 
section, "Subprograms") 

Exponentiation (**) 2nd 

Multiplication and Division (* and /) 3rd 

Addition and Subtraction (+ and -) 4th 


In addition, if two operators of the same hierarchy (with the 
exception of exponentiation) are used consecutively, the component 
operations of the expression are performed from left to right. 
Thus, the arithmetic expression A/B*C is evaluated as if the result 
of the division of A by B were multiplied by Cc. 


For example, the expression: 
(A*B/C*¥*#I+D) 
is effectively evaluated in the following order: 
a. A*B Call the result X (multiplication) (X/C**I+D) 
b. cC**I Call the result Y (exponentiation) (X/Y+tD) 


ce. X/Y Call the result Z (division) (Z+D) 
d. Z+D Final operation (addition) 


Note: This order of computation is used in determining the mode of 
an expression (see Table 2). 


For exponentiation the evaluation is from right to left. Thus, the 
expression: 


—A**B**C 


is evaluated as follows: 


a. B*¥*C Call the result Z 
b. A**Z Call the result Y 
Cs =x Final operation 


Use of Parentheses: Parentheses may be used in arithmetic expres- 


Sions, as in algebra, to specify the order in which the arithmetic 
operations are to be computed. Where parentheses are used, the 
expression within the parentheses is evaluated before the result is 
used. 


For example, the following expression: 


(B+ ( (A+B) *C) +A**2) 


is effectively evaluated in the following order: 


ae (A+B) Call the result X (B+ (X*¥C) +A**2) 
b. (X*C) Call the result Y (B+Y+A*¥*2) 

ce. BtyY Call the result W (W+A** 2) 

ad. A**2 Call the result Z (W+Z) 

e. WZ Final operation 
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7. j%Integer Division: When one integer is divided by another, the 
quotient is also an integer. If necessary, the result is truncat- 
ed. For example: 


5/2 


gives a quotient of 2. 


ARRAYS 


A FORTRAN array is a set of variables identified by a single variable 
name. A particular variable in the array may be referred to by its 
position in the array (e.g., first variable, third variable, seventh 
variable, etc.). Consider the array named NEXT which consists of five 
variables, each currently representing the following valves: 273, 41, 
8976, 59, and 2. 


NEXT(1) is the representation of 273 
NEXT(2) is the representation of 41 
NEXT(3) is the representation of &976 
NEXT(4) is the representation of 59 
NEXT(5) is the representation of 2 


Each variable in this array consists of the name of the array (i-e., 
NEXT) immediately followed by a number enclosed in varentheses, calied a 
subscript. The variables which comprise the array are called subscript- 
ed variables. Therefore, the subscripted variable NEXT(1) has the value 
273; the subscripted variable NEXT(2) has the value 41, etc. 


The subscripted variable NEXT(I) refers to the "“Ith" subscripted 
variable in the array, where I is an integer variable that may be 
asSSigned a value of 1, 2, 3, 4, or 5. 


To refer to the first element in an array, the array name must _ be 


subscripted. The array name itself does not represent the first 
element. 


Consider the following array named LIST consisting of two subscripts, 
the first ranging from 1 through 5, the second from 1 through 3: 


Columnl Column2 Column3 














Rowl 82 4 7 
Row2 12 13 14 
Row3 91 1 31 
Row4 24 16 10 
Row5 2 8 2 





Suppose it is desired to refer to the number in row 2, column 3; this 
would be: 


LIST (2,3) 


Thus, LIST (2,3) has the value 14 and LIST (4,1) has the value 24. 


18 


Ordinary mathematical notation might use LIST i,j to represent any 
element of the array LIST. In FORTRAN, this is written as LIST(I,J) 
where I equals 1,2,3,4, or 5 and J equals 1,2, or 3. 


As a further example, consider the array named COST consisting of 
three subscripts. This array might be used to store all the premiums 
for a life insurance applicant given (1) age, (2) sex, and (3) size of 
life insurance coverage desired. A code number could be developed for 
each statistic where IAGE represents age, ISEX represents sex, and ISIZE 
represents policy size desired. (See Table 4.) 


Table 4. Insurance Premium Codes 

(SS ee ee ere a ne ig a Pe eas 1 
| AGE | SEX | 
a aaa (insta Ss Ses SS ee { 
| | | 
| Age in Yrs. Code | Sex Code | 
| | | 
| | Male ISEX=1 | 
| d=: 5 IAGE=1 { Female ISEX=2 { 
| 6 - 10 IAGE=2 }—-—-—---—-~—-—-—-———-—-—-—---—--—--—-—-—--—-—--—--—--- { 
| 11 ~ 15 ITAGE=3 | POLICY SIZE | 
| 16 - 20 IAGE=4 -}--~-——-—-——--—~——~-~—~—-——--—-——-—-—~—~--—--—--——--- | 
| 21 25 IAGE=5 | Dollars Code | 
| 26 - 30 IAGE=6 | | 
| 31 35 IAGE=7 | 1,000 ISIZE=1 | 
i 36 4Q IAGE=8 | 2,000 ISIZE=2 | 
i 441 45 IAGE=9 | 3,000 ISIZE=3 | 
; 46 50 IAGE=10 | 5,000 ISIZE=4 | 
| P | 10,000 ISIZE=5 | 
| P | 25,000 ISIZE=6 | 
| : | 50,000 ISIZE=7 | 
| 96 100 IAGE=20 | 100,000 ISIZE=8 | 
ener ee ete Ey TTI ne eer 


Suppose an applicant were 14 years old, male, and desired a policy of 


$25,000. From Table 4, these statistics could be represented by the 
codes: 

IAGE= 3 (11 - 15 years old) 

ISEX=1 (male) 

ISIZE=6 ($25,000 policy) 


Thus, COST (3, 1, 6) represents the premium for a policy given the 
Statistics above. Note that "IAGE" can vary from 1 through 20, “ISEX" 
from 1 through 2, and “ISIZE" from 1 through 8. The number of 
subscripted variables in the array COST is the number of combinations 
that can be formed for different ages, sex, and policy Size available - 
a total of 20x2x8 or 320. Therefore, there may be up to 320 different 
premiums stored in the array named COST. 


The actual size (in storage locations) of the array COST depends upon 
the type of elements in the array. If each element in the array is 
real, the array size 1s 1280 storage locations; if each element is 
double precision, the array size is 2560 storage locations. 
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SUBSCRIPTS 


A subscript is a number used to refer to a particular variable within 
an array. A subscript may take one of several forms and there may be a 
maximum of three subscripts used with an array name. If more than one 
subscript is used, they must be separated by commas. All of the 
subscripts used with a particular array name must be enclosed in 
parentheses. The number of subscripts used must be equal to the number 
of dimensions in the array. 


A A AS NN NS SOR EY AOS AOE AES ANE AS ENN ED MD GND ED GRD ND SND GREE ONES SEU STD OR GEES AE GINS SED GED TEN ONY AUN COUNTY AR AED AED OA “NID AED ID ED SD ANY UD EID GENER AND SANIT OREN SN SENET REND AUNT SAND OD AERIS <ENUEE AEE <OUND AS SUNNY <I EE GID CEES SORTED 


| 

| Vv 

| ct 

i vtc' 

i v=cC" 

i Cc¥*V 

| c¥vtc' 
| c*v-c' 
| 

| 

| 

| 


Where: v represents an unsigned, nonsubscripted, integer variable. 


Whatever subscript form? is used, its evaluated result, as well as the 
intermediate result, must always be greater than 0 and less than or 
equal to 32,767. For example, when reference is made to the subscripted 
variable V(I-2), the value of I should be greater than 2 and less’ than 
or equal to 32,767. In any case, the evaluated result must be within 
the range of the array. 


Examples: 


Valid Subscripted Variables: 


ARRAY (IHOLD) 
NEXT (19) 
MATRIX (I-5) 
A(5*L) 
W(4*M+3) 


Invalid Subscripted Variables 


ARRAY (-T) (the subscript I may not be signed) 

COST (A+2) (A is not an integer variable unless defined as 
such by an Explicit specification statement) 

ARRAY (I+2.) (the constant within a subscript must be an 
integer) 

NEXT (-7*J) (the constant within a subscript must be unsigned) 

W(I(2)) (the subscript, I, may not be subscripted) 

LOT (0) (a subscript may never be nor assume a value of 
zero) 


1If more than one subscript form is used, the product of all subscripts 
must be less than or equal to 131,068 in Operating System FORTRAN IV (E) 
and less than or equal to 32,767 in the other three systems. 
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TEST (K*2) (if multiplication is indicated, the constant must 
precede the variable. Thus, TEST (2*K) is 
correct.) 

TOTAL (2+K) (if addition is indicated, the variable must pre- 
cede the constant. Thus, TOTAL (K+2) is correct.) 


DECLARING THE SIZE OF AN ARRAY 


The size of an array is determined by the number of subscripts of the 
array and the maximum value of each subscript. This information must be 
given for all arrays before using them in a FORTRAN program so that an 
appropriate amount of storage may be reserved. Declaration of this 
information is made by a DIMENSION statement, a COMMON statement, or by 
one of the Explicit specification statements; these statements are 
discussed in further detail in the section, “Specification Statements." 


ARRANGEMENT OF ARRAYS IN STORAGE 


Arrays are stored in ascending storage locations, with the value of 
the first of their subscripts increasing most rapidly and the value of 
the last increasing least rapidly. 


The array named A, consisting of one subscript which varies from 1 to 
5, appears in storage as follows: 


A(i) A(2) AC3) ACH) Al5) 
The array named B, consisting of two subscripts, whose first 


Subscript varies over the range from 1 to 5, and second varies from 1 to 
3, appears in ascending storage locations in the following order: 


B(1,1) B(2,1) B(3,1) BC4,1) B(5,1) 
Ba, 2) B(2,2) B(3,2) B4,2) Ess) 
B(1,3) B(2,3) B(3,3) B(4,3) B(5,3) 


Note that B(1,2) and B(1,3) follow in storage B(5,1) and B(5,2), 
respectively. 


The following list is the order of an array named C, consisting of 
three subscripts, whose first subscript varies from 1 to 3, second 
varies from 1 to 2, and third varies from 1 to 3: 


C4(1,.1, 1). CC2,172). C03, 1,40 CCl, 2,1) €(2,2,1) €G,23)) 





Cit ig 2) (CC2Z, be?) Clay ly. 2) Clig2s 2) C2722) C(3,242) 4 
CCL, 1,3): Cl2, 4,3) €Cl3,1,3) CO, 273): ClC2, 233) CO;2;3) 
Note that C(1i,1,2) and C(1,1,3) follow in storage C(3,2,1) and 
C(3,2,2), respectively. 
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ARITHMETIC STATEMENT 


a a aa a a a ec aa | 
| General Form | 
}----~~---------------------------------------------------------------- { 
| | 
J a=b | 
| , , | 
| Where: a is any subscripted or nonsubscripted variable. | 
| | 
i bis any arithmetic expression. | 
SE eee Ps ee Ue eure) psd wn SIP eg cern eo am OTN) ee Dey enOL Poh oe Oe PRE nO Ae cae Te Be RO mR RECs NE Peevey RoE J 


This FORTRAN statement closely resembles a conventional algebraic 
equation; however, the equal sign specifies replacement rather than 
equivalence. That is, the expression to the right of the equal sign is 
evaluated, and the resulting value replaces the current value of the 
variable to the left of the equal sign. 


Assume that the type of the following variables has been specified 
as: 


Variable Names Type 


I, J, W Integer variables 

A, B, D Real variables 

E Double-Precision variable 
F Real array 


Then the following examples illustrate valid arithmetic statements 
using constants, variables, and subscripted variables of different 
types: 


Statements Description 

A= B The value of A is replaced by the current value of B. 
W=B The value of B is truncated to an integer value, and 

this value replaces the value of W. 
A = iI The value of I is converted to a real value, and this 
: result replaces the value of A. 
I=I+#d The value of I is replaced by the value of I + 1. 
B = [*4*J+D The value of I is raised to the power J and result is 


converted to a real value to which the value of D is 
added. This result then replaces the value of B. 


A = BtD The most significant part of the product of B and D 

replaces the value of A. 
= I+E The value of I is converted to double precision and 

and added to E. The result of the addition is 
truncated from double precision to real and replaces 
the value of A. 

A = F(5,4) The value of F(5,4) replaces the value of A. 

E= 1 The value of I is converted to double precision, and 
this value replaces the value of E. 

J=E The value of E is truncated to an integer value, and 
this value replaces the value of J. 

E=A The value of A is converted to double precision, and 


this value replaces the value of E. 


CONTROL STATEMENTS 


Normally, FORTRAN statements are executed sequentially. That is, 
after one statement has been executed, the statement immediately 
following it is executed. This section discusses the statements that 
may be used to alter and control the normal sequence cf execution of 
Statements in the program. 


THE GO TO STATEMENTS 


The GO TO statements transfer control to the statement specified by 
number in the GO TO statement. Control may be transferred either 
unconditionally or conditionally. The GO TO statements are: 


1. The Unconditional GO TO Statement. 


2. The Computed GO TO Statement. 


Unconditional GO TO Statement 


GO TO xxXxxx 


Where: xXXxXxXX 1S an executable statement number. 


| | 
| | 
| | 
| | 
L J 


This GO TO statement causes control to be transferred to the 
statement specified by the statement number. Every subsequent execution 
of this GO TO statement results ina transfer to that same statement. 
Any executable statement immediately following this statement should 
have a statement number; otherwise, it can never be referred to or 
executed. 


Example: 


50 GO TO 25 
10 A=B+tC 
25 C = E**¥2 


@ 


Explanation: 


In the above example, every time statement 50 is executed, control is 
transferred to statement 25. | 
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Computed GO TO Statement 


i a 
fe ae Oe 


This statement causes control to be transferred to the statement 
numbered X44, X2se X3e+-+, OF Xn, Aepending on whether the current value 
of iis 1, 2, 3,..-, Or n, respectively. If the value of 1 is outside 
the allowable range, the next statement is executed. 


Example: 
GO TO (25, 10, 50, 7), ITEM 


50 A = BtC 


7 C = E**2+A 
25 L=C 


10 B = 21.3E02 


Explanation: 


In this example, if the value of the integer variable ITEM is 1, 
statement 25 will be executed next. If ITEM is equal to 2, statement 10 
is executed next, and so on. 


ADDITIONAL CONTROL STATEMENTS 


Arithmetic IF Statement 


De pe re ae eT ae eer Ot ee, ee ae ree 1 
| General Form | 
}---------------------------------------------------------------------- { 
| | | 
| IF (a) XasXaeX3 | 
| 
| Where: ais an arithmetic expression. | 
| | 
[ XaeXaeX3 are executable statement numbers. | 
a i a a ee ll sc Sa hc me mes eae a i a a a ca le Jj 


This statement causes control to be transferred to the statement 
numbered xX4,X2, OF X3 when the value of the arithmetic expression (a) is 
less than, equal to, or greater than zero, respectively. The first 
executable statement following the arithmetic IF statement should have a 
Statement number; otherwise, it can never be referred to or executed. 


Exanpie: 


IF (A(J,K)**3-B)10, 4, 30 


II 
) 
* 
# 
N 


30 C 


10 E (F*B)/D+1 


Explanation: 


In the above example, if the value of the expression (A(J,K)¥**3-B) is 
negative, the statement numbered 10 is executed next. If the value of 
the expression is zero, the statement numbered 4 is executed next. If 
the value of the expression is positive, the statement numbered 30 is 
executed next. 


DO Statement 


re a a eg ee ee gg ee ee ee gee a a ee Pe ee EN yee ge 1 
| General Form | 
fm a a i a ee | 
| End of DO Initial Test ( 
| Range Variable Value Value Increment | 
| DO x 2; = Ma, Ma, M3 
| | 
| Where: x iS an executable statement number that is not defined | 
| before the DO statement. i 
| | 
{ iis a nonsubscripted integer variable. | 
| M4, Ma, M3, are either unsigned integer constants greater | 
{ than zero or unsigned nonsubscripted integer variables whose | 
i value is greater than zero. Mz May not exceed 234-2 in | 
i value. m3, is optional; if it is omitted, its value is | 
i assumed to be 1. In this case, the preceding comma must | 
| also be omitted. | 
eee eee RENE Noa Se RRO OER TON Soe A SO ACTEM Me aS SNORE EOE ON MMCTE Nate nN ene Ree ee eR Tes Se ee Te ee ee eT 


The DO Statement is a command to execute repeatedly the statements 
that follow, up to and including the statement numbered x. The first 
time the statements in the range of the DO are executed, i is 
initialized to the value m,; each succeeding time i is increased by the 
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value m3- When, at the end of the iteration, i is equal to the highest 
value that does not exceed mz, control passes to the statement following 
the statement numbered x. Thus, the number of times the statements in 
the range of the DO is executed is given by the expression: 


where the brackets represent the largest integral value not exceeding 
the value of the expression. If mz is less than m, the statements in 
the range of the DO are executed once. Upon completion of the DO, the 
DO variable is undefined, and should not be used until redefined (e.g., 
in a READ list). 


There are several ways in which looping (repetitively executing the 
Same statements) may be accomplished when using the FORTRAN language. 


Example 1: 


Assume that a manufacturer carries 1,000 different machine parts in 
Stock. Periodically, he may find it necessary to compute the amount of 
each different part presently available. This amount may be calculated 
by subtracting the number of each item used, OUT(I), from the previous 
Stock on hand, STOCK(I). 

5 I=0 

10 I=I+1 

25 STOCK (I)=STOCK(I)- OUT(T) 

15 IF (I-1000) 10,30,30 

30 A=Btc 


Explanation: 


The three statements (5, 10, and 15) required to control the 
previously shown loop could be replaced by a single DO statement as 
follows: 


DO 25 I = 1,1000 
25 STOCK(I) = STOCK (I)-OUT(T) 
30 A = BtC 


In the above code, the DO variable, I, is set to the initial value of 1. 
Before the second execution of statement 25, I is increased by the 
increment, 1, and statement 25 is again executed. After 1000 executions 
of the DO loop, I equals 1000. Since I is now equal to the highest 
value that does not exceed the test value, 1000, control passes out of 
the DO loop and statement 30 is executed next. Note that the DO 
variable I is now undefined; its value is not necessarily 1000 or 1001. 
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Example 2: 


DO 25 I=1, 10, 2 
15 J = I+K 
25 ARRAY(J) = BRAY (J) 
30 A=B+#t#C 


Explanation: 


In the preceding example, statement 25 is the end of the range of the 
DO loop. The DO variable, I, is set to the initial value of 1. Before 
the second execution of the DO loop, I is increased by the increment, 2, 
and statements 15 and 25 are executed a second time. After the fifth 
execution of the DO loop, I equals 9. Since I iS now egual to the 
highest value that does not exceed the test value, 10, control passes 
out of the DO loop and statement 30 is executed next. Note that the DO 
variable I is now undefined; its value is not necessarily 9 or 11. 


Programming Considerations in Using a DO Loop 


1. The indexing parameters of a DO statement (i, Mm, Mig, M3) may not 
be changed by a statement within the range of the DO loop. 


2. There may be other DO statements within the range of a DO 
statement. All statements in the range of the inner DO must be in 


the range of the outer DO. A set of DO statements satisfying this 
rule is called a nest of DO's. 


Example 1: 
DO 50 T= 1, 4 


A(I) = B(I) **2 


DO 50 J=1, 5 Range of 
Range of Outer DO 
50 C(J+1) = A(T) Inner DO 
Example 2: 


DO 10 INDEX = L, M 
N = INDEX + K 
DO 15 J = 1, 100, 2 Range of 
Range of Outer DO 
15 TABLE(J) = SUM(J,N)-1 Inner DO 
10 B(N) = A(N) 
3. A transfer out of the range of any DO loop is permissible at any 


time; a transfer into the range of a DO loop is permissible only as 
described in item 4. 
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4, When a transfer is made out of the range of an innermost DO _ loop, 
transfer back into the range of that DO loop is allowed if and only 
if none of the indexing parameters (i,M,,M2,M3) are changed outside 
the range of the DO 


Example: 
DO 


207 





Explanation: 


In the preceding example, the transfers specified by the numbers 
1,2, and 3 are permissible, whereas those specified by 4,5, and 6 
are not. 


5- The indexing parameters (1,%,,M2,m3) may be changed by statements 
outside the range of the DO statement only if no transfer is made 
back into the range of the DO statement that uses those parameters. 


6. The last statement in the range of a DO loop (statement x) must be 
an executable statement, not of the form GO TO, PAUSE, STOP, 
RETURN, Arithmetic IF, or another DO. 


7. The use of, and return from, a subprogram from within any DO loop 
in a nest of DOs is permitted. 


CONTINUE Statement 


| General Form | 


CONTINUE is a dummy statement which may be placed anywhere in the 
source program without affecting the sequence of execution. It may be 
used as the last statement in the range of a DO in order to avoid ending 
the DO loop with a GO TO, PAUSE, STOP, RETURN, Arithmetic IF or another 
DO statement. 
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Example 1: 


DO 30 I = 1, 20 
7 IF (A(I)-B(I)) 5,30,30 
5 ACI) =A(I) +1.0 

B(I) = B(I) -2.0 


@ 


GO TO 7 
30 CONTINUE 
40 C = A(3) + B(7) 


Explanation: 


In the preceding example, the CONTINUE statement is used as the last 
Statement in the range of the DO in order to avoid ending the DO loop 
with the statement GO TO 7. 


Example 2: 


DO 30 I=1, 20 

IF (A(I)-B(I))5, 40,40 
5 A(I) = C(I) 

GO TO 30 
40 A(I) = 0.0 
30 CONTINUE 


Explanation: 


In Example 2, the CONTINUE statement provides a branch point enabling 
the programmer to bypass the execution of statement 40. 


PAUSE Statement 


EE AE rE Ee TS ORL ERS EE A NRE ET UE CINDER OD EEN EE NET A SEEN NED OEE AD aD ED ER GOS A ED OD CD CUNY RE ART OD ERE LED ERY GED GRE ETD AES GEN WORD UAL CD SIMA SETS SEEN RED EERE ENN ENED SND IONE SEED GENE CSD GREED UTD EE GARNI tee celeEED GG GREED eee ae ED 


| General Form | 
|---------------------------------------------------------------------- { 
| PAUSE | 
| PAUSE n | 
| 
| Where: n is an unsigned 1 through 5 digit integer constant. [ 
Ba ee een ee eee a 4 


Information is displayed and the program waits until operator 
intervention causes it to resume execution, starting with the next 
Statement after the PAUSE statement. The particular form of the PAUSE 
statement used determines the nature of the information that is 
displayed. The PAUSE statement causes PAUSE 00000 to be displayed. If 
n is specified, PAUSE n is displayed. 
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STOP Statement 


| General Form | 


This statement terminates the execution of the object program and 
displays n if specified. 


END Statement 


| END | 


Ve oe AD A ETD NS SE AS ED ES TS ET ND ERTS EAD SEED ERS ED ND GENS GEREED ED REND EEN EEDTD CANNED SERED QURTS EES SERED AED EENED CETUS GERESD GREED EEE GRE GERD IT END TD EEE GEES CED GERD OUT qu Gee <nteD ame CREE <u qae dee eube ner CRE: | 


The END statement is a nonexecutable statement that defines the end 
of a source program or source subprogram for the compiler. Physically, 
it must be the last statement of each program or subprogram, and it may 
not be continued. 
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INPUT/OUTPUT STATEMENTS 


The input/output statements enable a user to transfer data, belonging 
to a named collection of data, between input/output devices (such as 
disk units, card readers, and magnetic tape units) and internal storage. 
The named collection of data is called a data set and is a continuous 
string of data that may be divided into FORTRAN records. 


A data set is referred to by an integer constant or integer variable. 
Formerly, this reference was called a symbolic unit number. However, 
because the reference is to the data rather than any specific device, 
this number is called the data set reference number. 


There are two types of I/0 statements: sequential I/O statements 
(available in all Basic FORTRAN IV systems) and direct access I/0 
statements (not available in Basic Programming Support FORTRAN IV). The 
sequential statements provide facilities for the sequential selection 
and placement of data. These statements are device independent because 
a given statement may be applicable to a data set On any number of 
devices or device types. 


The direct access I/O statements provide facilities for the selection 
and placement of data in an order specified by the user. These 
Statements are only valid when the data set will be or is already 
resident on a direct access storage device. 


SEQUENTIAL INPUT/OUTPUT STATEMENTS 


There are five sequential I/O statements: READ, WRITE, END FILE, 
REWIND, and BACKSPACE. The READ and WRITE statements cause transfer of 
records of sequential data sets. The END FILE statement defines the end 
of a data set; the REWIND and BACKSPACE statements control the 
positioning of data sets. 


In addition to these five statements, the FORMAT statement, although 
it is not an I/O statement, is used with certain forms of the READ and 
WRITE statements. The FORMAT statement specifies the form in which the 
data is to be transmitted and allows the user to divide a data set into 
FORTRAN records. 


Even though the I/O statements are device independent, the original 
source or the ultimate destination of the data being transferred 
influences the specification of the records and data formats. There- 
fore, subsequent examples are in terms of card input and print-line 
output unless otherwise noted. 
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READ STATEMENT 


ee alg cee ee ee are ee ne paw ee SS pee Sey pe erg ee ge ee eee ae eae 


{| General Form | | | 


READ(a, b) list 


Where: ais an unsigned integer constant or an integer variable 
that represents a data set reference number. 


describes the data being read. 


list is a series of variable or array names, which may be 
indexed and must be separated by commas. These names 
specify the number of items to be read and the locations in 
storage into which the data is placed. 


CLD) OLED SOY RES EDN ETT SOLED AEE IE ATED CD ED INT LE EE SE UP SY ED EET EGE RED ORR NED COUR GREE GE EE SE AEEED CER CRED CES eRe cae Re EE oD SEE EE ER eo GS oP <a were: Es <e eD ee eee ee -ENer quEse eines emmy etre ere come «eee eu ane Gee wR eer anit ee ieee one nD 


i | 
{ | 
i | 
| | 
| | 
| ! 
| b is the statement number of the FORMAT statement that | 
| | 
| 
| | 
| 
| | 
| | 
L J 


The READ statement may take many forms. Either the list parameter or 
the b parameter may be omitted. 


The basic forms of the READ statement involve formatted and unformat- 
ted data. They are respectively: 


READ (a, b) list 
READ (a) list 





The Form READ (a,b) list 


This form is used to read data from the data set associated with a 
into the locations in storage specified by the variable names in the 
list. The list, used in conjunction with the specified FORMAT statement 
b (see the section, “FORMAT statement"), determines the number of items 
(data) to be read, the locations, and the form the data will take in 
storage. 


Example 1: 


Assume that the variables A, B, and C have been declared as integer 
variables. 


75 FORMAT (110, I8, I9) 


READ (J, 75) A, B, C 
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Explanation: 


The above READ statement causes input data from the data set 
associated with data set reference number J to be read into the 
locations A, B, and C as specified by FORMAT statement 75. That is, the 
first 10 positions of the record are read into storage location A; the 
next 8 positions are read into storage location B; and the next 9 
positions are read into storage location C. 


The list may be omitted from the READ (a,b)list statement. In this 
case, a record is skipped or data is read from the data set associated 
with a into the locations in storage occupied by the FORMAT statement 
numbered b. 





Example 2: 


98 FORMAT ("HEADING") 


Explanation: 


The above statements would cause the characters H, E, A, D, I, N, and 
G in storage to be replaced by the next 7 characters in the data _ set 
associated with data set reference number 5. 


Example 3: 


98 FORMAT (110, "HEADING') 


READ (5,98) 


Explanation: 


The above statements would cause the next record in the data set 
associated with data set reference number 5 to be skipped. No data is 
transferred into internal storage because there is no list item which 
corresponds with format code I[10. 


The Form READ (a) list 


The form READ (a) list of the READ statement causes binary data 
(internal form) to be read from the data set associated with a into the 
locations of storage specified by the variable names in the list. Since 
the input data is always in internal form, a FORMAT statement is not 
required. This statement is used to retrieve the data written by a 
WRITE (a) list statement. 
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Example 1: 


READ (5) A, B, C 


Explanation: 


This statement causes the binary data from the data set associated 
with data set reference number 5 to be read into the storage locations 
specified by the variable names A, B, and C. 


The list may be omitted from the READ (a) list statement. In this 
case, a record is skipped. 


Example 2: 
READ (5) 


Explanation: 


The above statement would cause the next record in the data set 
associated with data set reference number 5 to be skipped. No data is 
transferred into internal storage. 


Indexing I/O Lists 


Variables within an I/O list may be indexed and incremented in the 
Same manner as those within a DO statement. These variables and their 
indexes must be included in parentheses. For example, suppose it is 
desired to read data into the first five positions of the array A. This 
may be accomplished by using an indexed list as follows: 


15 FORMAT (F10.3) 


READ (2,15) (A(I), I=1,5) 
This is equivalent to: 


15 FORMAT (F10.3) 


DO 12 I =1,5 
12 READ (2,15) A(T) 


As with DO statements, a third indexing parameter may be used to 
specify the amount by which the index is to be incremented at each 
iteration. Thus, 

READ (2,15) (A(T), J=1,10,2) 
causes transmission of values for A(1), A(3), A(5), A(7), and A(9). 
Furthermore, this notation may be nested. For example, the state- 


ment: 


READ (2,15) ((C(I,d),D(1I,d),J=1, 3), I=1, 4) 
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would transmit data in the following order: 


C(L,y1), D,1), CU, 2)4 DOy2), CO, 6 DO, 3),% 
C(2,1), D(2,1), C(2,2), D(2,2), C(2,3), D(2,3), 
C(3,1), D(3,1), C(3,2), D(3,2), C(3,3), DG,3), 
C(4,1), DC4,1), C(4,2), DC4,2), C(4,3), DC4,3). 


Since J is the innermost index, it varies more rapidly than I. 
As another example, consider the following: 


READ (2,25) I, (C(J),J=1,1) 


The variable I is read first and its value then serves aS an index to 
specify the number of data items to be read into the array C. 


If it is desired to read data into an entire array, it is not 
necessary to index that array in the I/O list. For example, assume that 
the array A consists of one subscript ranging from ito 10. Then the 
following READ statement referring to FORMAT statement numbered 5: 


READ (2,5) A 
would cause data to be read into A(1), A(2),...,A(10). 


The indexing of I/0 lists applies to WRITE lists as well as READ 
lists. 


WRITE STATEMENT 


| General Form 


atlas AP RIED AD GORA ERI aOR! RMI ARIES ete CERES EERE aii AERTS ONE EN ER ER OOO cot eR END CED OCD MEDRN SED aEE AED “RD deme SEH we ES aoe GRD GEOR coe GSS aE eemdy emran tle came ene wD ED Ee OME Ee wo comme eee ee ee ee ene ID woo mee ee ome eA ean cee eee coe eS CE eo 


WRITE (a, b) list 


Where: ais an unsigned integer constant or an integer variable 
that represents a data set reference number. 


b is the statement number of the FORMAT statement that 
describes the data being written. 


list is a series of variable or array names, which may be 
indexed and must be separated by commas. These names 
specify the number of items to be written and the locations 
in storage from which the data is taken. 


QPP sone cone eee eS se ee CED semen ED GE gee Gee ee 


The WRITE statement may take many different forms. For example, the 
list or the parameter b may be omitted. 


The basic forms of the WRITE statement involve formatted and 
unformatted data. They are: 
WRITE (a,b) list 
WRITE (a) list 
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The Form WRITE (a,b) list 


This form is used to write data in the data set associated with a 
from the locations in storage specified by the variable names in the 
list. The list, used in conjunction with the specified FORMAT statement 
b, determines the number of items (data) to be written, the locations, 
and the form the data will take in the data set. 


Example 1: 


In the following example, assume that the variables A, B, and C have 
been declared as integer variables. 


75 FORMAT (110, 18, 19) 


1@ 


WRITE (J, 75) A, B, C 


Explanation: 


The above WRITE statement causes output data to be written in the 
data set associated with the data set reference number J, from the 
locations A, B, C, as specified by FORMAT statement 75. That is, the 10 
rightmost digits in A are written in the data set associated with the 
data set reference number J; the next 8 positions in the data set will 
contain the 8 rightmost digits in B; and the next 9 positions in the 
data set will contain the 9 rightmost digits in Cc. 


The list may be omitted from the WRITE (a,b) list statement. In this 
case, a blank record is inserted or data is written in the data set 
associated with a from the locations in storage occupied by the FORMAT 
statement b. 


Example 2: 
98 FORMAT (* HEADINGS ) 


WRITE (5,98) 


The above statements would cause a blank and the characters H, EF, A, 
D, I, N, and G in storage to be written in the data set associated with 
data set reference number 5. 


Example 3: 
98 FORMAT (110, ‘HEADING') 


WRITE (5,98) 


Explanation: 


The above statements would cause a blank record to be written in the 
data set associated with data set reference number 5. No data is 
transferred into the data set because there is no list item which 
corresponds with the format code I10. | 


36 


The Form WRITE (a) list 


The WRITE (a) list form of the WRITE statement causes binary data 
(internal form) from the locations of storage specified by the variable 
names in the list to be written in the data set associated with a. 
Since the output data is always in internal form, a FORMAT statement is 
not required. The READ (a) list statement is used to retrieve the data 
written by a WRITE (a) list statement. 





Example: 
WRITE (5) A, B, C 


Explanation: 


The statement causes the binary data from the locations specified by 
the variable names A, B, and C to be written in the data set associated 
with data set reference number 5. | 


FORMAT STATEMENT 
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| General Form 


XXXXX FORMAT (C4 ,Ca,z eee eCn/C,' peas eee *Cn'/. « -) 





Day 
> 
@ 
K 
1) 


XXXxXx is a statement number (1 through 5 digits). 


CasCaneeeeCn and Cy", Ca',---,Cyn*' are format codes which may 
be delimited by one of the separators: comma, slash, or 
parenthesis. These codes specify the length, decimai point 
(if any), and position of the data in the data set. 


/ may be used to separate FORTRAN records. 


ff me ee ee ee ee ee eee ee eee ee oe oe 


The FORMAT statement is used in conjunction with the READ and WRITE 
statements in order to specify the desired form of the data to be 
transmitted. The form of the data is varied by the use of different 
format codes. 


The format codes are: 


I - to transfer integer data 

F - to transfer real or double precision data that does not contain 
a decimal exponent 

E or D - to transfer real or double precision data that contains an 
E or D decimal exponent, respectively 

A - to transfer character data 

Literal - to transfer a string of alphameric and special characters 

H - to transfer literal data 

~ to either skip data when reading or insert blanks when writing 

- to specify the position in a FORTRAN record where transfer of 

data is to start 

to specify a scale factor 


i * eB 
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Any number used in a FORMAT statement, except the statement number, 
scale factor, or any literal, must be less than or equal to 255. The 
scale factor must be less than or equal to 127. 


USE OF THE FORMAT STATEMENT: This section contains general information 
on the FORMAT statement. The points discussed below are illustrated by 
the examples that follow. 


1. FORMAT statements are nonexecutable and may be placed anywhere in 
the source program. 


2. A FORMAT statement may be used to define a FORTRAN record as 
follows: 


a. If no slashes or additional parentheses appear within a FORMAT 
statement, a FORTRAN record is defined by the beginning of the 
FORMAT statement (left parenthesis) to the end of the FORMAT 
statement (right parenthesis). Thus, a new record is read when 
the format control is initiated (left parenthesis); a new 
record is written when the format control is terminated (right 
parenthesis). 


Example: 


XXXXX FORMAT (----, ----, ---~-) 


~--corresponds to 1 
FORTRAN record 


b. If slashes appear within a FORMAT statement, FORTRAN records 
are defined by the beginning of the FORMAT statement to the 
first slash in the FORMAT statement, from one slash to the next 
succeeding slash, or from the last slash to the end of the 
FORMAT statement. Thus, a new record is read when the format 
control is initiated, and thereafter a record is read upon 
encountering a slash; a new record is written upon encountering 
a Slash or when format control is terminated. 


Example: 
XXxXxXxX FORMAT (----/ <----/ ----) 


<---> <---> <---> 


a ee a ee each corresponds to 
1 FORTRAN record 


c. If more than one level of parentheses appear within a FORMAT 
statement, a FORTRAN record is defined by the beginning of the 
FORMAT statement to the end of the FORMAT statement. At this 
point, the definition of the FORTRAN record continues at the 
first-level left parenthesis that precedes the end of the 
FORMAT statement. | 
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Example: 


0 al 1 1 1 0 
XXX. FORMAT Cee]: (+]-) S25 (S2=)° ==) 
aim a aaa as a > 
| 
| <-=s->-= > 
| | 
| 
| | 
t—___———___—~—~~—~—~——~——~— each corresponds to 


1 FORTRAN record 


When defining a FORTRAN record by a FORMAT statement it is 
important to consider the original source (input) or ultimate 
destination (output) of the record. For example, if a FORTRAN 
record is to be punched for output, the record should not be 
greater than 80 characters. For input, the FORMAT statement should 
not define a FORTRAN record longer+ than the record referred to in 
the data set. 


3. Blank output records may be introduced or input records may be 
skipped by using consecutive slashes (/) in a FORMAT statement. If 
there are n consecutive slashes at the beginning or end of a FORMAT 
statement, n input records are skipped or n blank records are 
inserted between output records, respectively. If n consecutive 
Slashes appear anywhere else in a FORMAT statement, the number of 
records skipped or blank records inserted is n-1. For example, the 
Statements: 


10 FORMAT (///1I6) 


@ 


READ (INPUT,10) MULT 


cause three records to be skipped on the data set associated with 
INPUT before data is read into MULT. 


The statements, where "x" is a carriage control character (see, 
"Carriage Control"): 


15 FORMAT (*x',15,////'x',F5.2,12//) 


WRITE (I0UT,15) K,A,J 
result in the following output: 


integer 
(blank line) 
(blank line) 
(blank line) 
Real, Integer 
(blank line) 
(blank line) 


1In Basic Programming Support FORTRAN IV, Disk Operating System FORTRAN 
IV, and Tape Operating System FORTRAN IV, a maximum of 255 characters 
per record may be transmitted. In Operating System FORTRAN IV (E), the 
maximum depends upon the device (see the FORTRAN IV (E) Programmer's 
Guide). 
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Successive items in an I/O list are transmitted according to 
successive format codes in the FORMAT statement, until all items in 
the list are transmitted. If there are more items in the list than 
there are codes in the FORMAT statement, control transfers to the 
preceding left parenthesis of the FORMAT statement and the same 
format codes are used again with the next record. If there are 
fewer items in the list, the remaining format codes are not used. 
For example, suppose the following statements are written in a 
program: 


10 FORMAT (F10.3,8£12.4,F12.2) 


WRITE (3,10) A,B,C,D,E,F,G 


The following table shows the data transmitted in the column on the 
left and its corresponding format code. 


Data Transmitted Format Codes 

A Fi0.3 first data 

B Ei2.4 record 

Cc F12.2 

D F10.3 second data 

E E12.4 record 

F F1i2.2 

G F10.3 third data 
record 


A format code may be repeated as many times as desired by preceding 
the format code with an unsigned integer constant. Thus, 


(2F10.4) 
is equivalent to: 
(F10.4,F10.4) 


A limited one-level parenthetical expression is permitted to enable 
repetition of data fields according to certain format codes within 
a longer FORMAT statement. If a multiline listing is desired such 
that the first two lines are to be printed according to a special 
format and all remaining lines according to another format, the 
last format code in the statement should be enclosed in a second 
pair of parentheses. For example, in the statement: 


FORMAT (‘'x',12,F3.1/'x',F10.8/('x'" ,3F5.1)) 


If more data items are to be transmitted after the format codes 
have been completely used, the format repeats from the last left 
parentheses. Thus, the printed output would take the following 
form: 


F10.8 


FP5.1,F5.1,F5.1 
F5<1,F5.1,F5.1 


As another example, consider the following statement: 


FORMAT (*x',1I2/2('x',13,F6.1) ,F9.7) 


If there are 13 data items to be transmitted, then the printed 
output on a WRITE statement would take the following form: 


I2 

3, F6614" X47 13 7, Foci ro. 7 
13,F6.1,;°x*,13,F6.1,F9.7 
I3,F6.1 


7. When transferring data on input or output, the type of format code 
used, type of data, and type of variables in the I/O list should 
correspond. 


8. In the following examples, the output is shown as a printed line. 
A carriage control character ‘x', (see, "Carriage Control") is 
Specified in the FORMAT statement but does not appear inthe first 
print position of the print line. This carriage control character 
appears as the first character of the output record on any I/0 
medium other than the printed line. 


Numeric Format Codes (I,F,E,D) 


Four types of format codes are available for the transfer of numeric 
data. These are specified in the following form: 


ih Saree ie a i Ee Se 
| General Form 

alw 

aFw.d 

aEw.d 

aDw.d 

Where: ais optional and is an unsigned integer constant used _ to 


denote the number of times the same format code is repeti- 
tively referenced. 


I,F,E,D are format codes. 


w is an unsigned integer constant that is less than or equal 
to 255 and specifies the number of characters of data. 


d is an unsigned integer constant that specifies the number 
of decimal places to the right of the decimal point, i.e., 
the fractional portion. 


ROO RG TS ERED ERED EES COTS EER ENS ED EE EE RD NED NE ED GOED ED EN TED ERED ANU RD GUD ENE GENES GOT EN ERED RD AOR SURED NED SEDONA SE AD AEE TD EES GED SED QE ED GET ETE GED SEUSS ETD ATED GEEKS CORED ULE CRED ENS GES COED OED GED eee cee, GEUEP NED RED ED ERD cee RE ED ee Ce 


Pr — oe eee oe cee ee eee ee ee eee ee ee ee ee ee ee ee ee 


For purposes of simplification, the following discussion of format 
codes deals with the printed line. The concepts developed apply to all 
input/output media. 


I Format Code 


The I format code is used to transmit integer data. If the quantity 
is negative, the position preceding the leftmost digit contains a minus 
Sign. In this case, an additional position should be specified in w for 
the minus sign. If the number of characters is less than w, on input, 
leading blanks are not’ significant, embedded and trailing blanks are 
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treated as zeros; on output, the leftmost print positions are filled 
with blanks. If the number of characters to be transmitted is greater 
than w, on output, asterisks are given. The magnitude of the data to be 
transmitted must not exceed the maximum magnitude for an integer 
constant. | 


The following examples show the internal value for each of the 
quantities on the left, according to format code I3: (b represents a 
blank): | 


External Value Internal Value 
345 345 
bb4 004 
Wb3 403 
43b 430 


The following examples show how each of the quantities on the left is 
printed according to the format code I3: 


Internal Value Printed Value 


721 721 

-721 721 (incorrect because of insufficient 
specification) 

-12 -12 

568114 se (incorrect because of insufficient 
specification) 

0 bb0 

-5 b-5 

9 bb9 


F Format Code 


The F format code is used in conjunction with the transferral of real 
or double precision data that does not contain a decimal exponent. For 
F format codes, w is the total field length reserved and d is the number 
of places to the right of the decimal point (the fractional portion). 
The total field length reserved must include sufficient positions for a 
Sign (if any) and a decimal point. The sign, if negative, is printed. 
The magnitude of the data to be transmitted must not exceed the maximum 
magnitude for a real or double precision constant. 


The integer portion of the number is handled in the same fashion as 
numbers transmitted by the I format code. If excessive positions are 
reserved by d, zeros are filled in on the right. If insufficient 
positions are reserved by d, on input the fractional portion is 
truncated from the right and on output the fractional portion is 
rounded; truncating and rounding occur at the dth position. 


The following examples show the internal value for each of the 
quantities on the left, according to format code F5.2: 


External Value Internal Value 
32.46 32.46 
-8.4 -8.40 
13.568 13.57 


The following examples show how each of the quantities on the left is 
printed according to the format code F5.2: 
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Internal Value Printed Value 


12.17 R247 

—-41.16 41.16 (incorrect because of insufficient 
specification) 

-.2 -0.20 

7.3582 b7.36 

~1. ~1.00 

9.03 b9. 03 

187.64 * &  & (incorrect because of insufficient 
specification) 


E and D Format Codes 


The E and D format codes are used in conjunction with the transferral 
Of real or double precision data that contains an E or D decimal 
exponent, respectively. For E and D format codes, the fractional 
portion is again indicated by d. 


The w includes field d, spaces for a Sign, the decimal point, plus 
four spaces for the exponent. For output, Space for at least one digit 
preceding the decimal point should be reserved. In general, w should be 
at least equal to dt?. 


The exponent is the power of 10 by which the number must _ be 
multiplied to obtain its true value. The exponent is written with a D 
or an E, followed by a space for the sign and two spaces for the 
exponent (maximum value is 75). 


The following examples show the internal value for each of the 
quantities on the left, according to the format codes (D10.3/E10.3) 


External Value Internal Value 
02.380+02 238. 
-0.001E+03 = 

-7.654D-06 ~.000007654 
1.EF2 100. 

4.673 “4.673 


The following examples show how each of the quantities on the left is 
printed, according to the format codes (D10.3/E10.3): 


Internal Value Printed Value 
238. bO. 238Db03 
-.002 -0.200E-02 
-00000000004 bO.400D-10 
-21.0057 -Q0.210Eb02 (Last three digits of accuracy 
lost because of insufficient 
specification) 


When reading input data, the start of the exponent field must either 
be marked by an E or D, or, if that is omitted, by a + or - sign (not a 
blank). Thus, E2,E+2,+2,+02,E02, and E+02 all have the same effect and 
are permissible decimal exponents for input. 


Numbers for E, D, and F format codes need not have their decimal 
point punched. If it is not present, the decimal point is supplied by 
the d portion of the format code. If it is present in the card, its 
position overrides the position indicated by the d portion of the format 
code. 
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A Format Code 


SE ND ENED AE A EE TE EE ED NS ET ES GENET SET EN ETD TS ENT NED GAT ENE CANS ERED SEED ALND ED ERA AED A SED GN SILI ND SENS ID ID AN ED AD NG A ES ED EOS IE SENS ERED GTR ENP REED GENTE ERED OYRED GEED NES ETD ERED SEED CEE GERD GED GED GUD SHEED CIES eoTED GIUED <OD aE ene 


| General Form 


SOUR AIEEE ARUN IER EN CT AER IEE OEE CINE OR TIES CREEL CE GURL AER LEY IND VAC ERED SIGNS ARERS CEA IED ED AEE OED CHEERS A RD OOS ETM GR A GERD sO ET EEE ES aE EE cite clutenD siamo OSHS Ee CER oie “eRe <PaE ates GE eee my te tee ee SO EY GRD eS Ee GD ele eee aoe VE cme aD ee 


is 
< 


Where: ais optional and is an unsigned integer constant used to 
denote the number of times the same format code is repeti- 
tively referenced. 


w is an unsigned integer constant that is less than or equal 
to 255 and specifies the number of characters of data, 
including blanks. 


ASD ARES FREES NOS ERIN SPEER: GENES GREED GERD SOE URE CET GOES QE) ED EAD “GEO GED EES ERE DE AS OD ES GERD ERED GEE GUND Ga Ge a GE ee ee ee cee eee ese <r <CeD e atee GS Ee ei ee ee ee ee ee ee ee ee oe oe ee oo ee ee ee 


QQ ee ee ee ee ee a a a a 


The format code Aw is used to read or write character data. If w is 
equal to the number of characters corresponding to the length of each 
item in the I/O list, w characters are read or written. 


On input, if wis less than the length of the storage reserved for 
each item in the I/O list, w characters are read and the remaining 
Yight-most characters in the item are replaced with blanks. If wis 
greater than the length, the number of characters equal to the 
difference between w and the length are skipped, and the remaining 
characters are read. 


On output, if w is less than the length of the storage reserved for 
each item, the printed line will consist of the left-most w characters 
of the item. If wis greater than the length the printed line will 
consist of the characters right-justified in the field and be preceded 
by blanks. Therefore it is important to always allocate enough storage 
to handle the characters being written (see the section "Specification 
Statements"). 


Example 1: 


Assume that the array ALPHA consists of one subscript ranging from 1 
through 20. The following statements could be written to “copy” a 
record from one data set to another whose ultimate destination is a card 
punch. 


10 FORMAT (20A4) 


READ (5,10) (ALPHA(I),I=1,20) 


WRITE (6,10) (ALPHA(I) ,I=1,20) 


G4 


Explanation: 


In this example, the REFAD statement would cause 20 groups of 
characters to be read from the data set associated with data set 
reference number 5. Each group of four characters would be placed into 
one of the 20 positions in storage starting with ALPHA(1) and ending 
With ALPHA(20). The WRITE statement would cause the 20 groups of four 
characters to be written on the data set associated with data set 
reference number 6. 


Example 2: 


As another example, consider all the variable names in the list of 
the following READ statement to have been explicitly specified as REAL 
and the array CONST to have been specified as having one subscript 
ranging from 1 through 10. Then assuming the following input data is 
associated with data set reference number 5, 


ABCDE.. .XY¥Z$1234567890b 


where ... represents the alphabetic characters F through W and b means 
a blank, the following statements could be written: 


10 FORMAT (27A1,10A1,A1) 
20 FORMAT ('x',6(7A1,5X)) 


e 


READ (5,10) 


A,B,C,D,E,F,G,H,I, 
1 J,K,L,M,N,0,P,9,R, 
2 ST, U,V iW,X,¥,2,5, 
3 (CONST (IND), IND=1, 10), BLANK 
DO 50 INDEX = 1,5 
WRITE (6, 20)G,R,0O,U,P,BLANK,CONST(INDEX), 
1 B, L,O,C,K, BLANK, CONST (INDEX), 
2 F,1,E,L,D, BLANK, CONST (INDEX), 
3 G,R, O,U, P, BLANK, CONST (INDEX+5), 
4 B,L,O,C,K, BLANK, CONST (INDEX+5), 
5 F, 1,E, L,D, BLANK, CONST (INDEX+5) 


a 


50 CONTINUE 


Explanation: 


The READ statement would cause the 37 alphameric characters and the 
blank in the data set associated with data set reference number 5 to be 
placed into the storage locations specified by the variable names in the 
READ list. Thus, the variables A through Z receive the values A through 
Ze, respectively; the variable $ receives the value $; the numbers 1 
through 9, and 0, are placed in the ten fields in storage starting with 
CONST(1) and ending with CONST(10); and the variable BLANK receives a 
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blank. The WRITE statement within the DO loop would cause the following 
heading to be printed. A subsequent WRITE statement within the DO loop 
could then be written to print the corresponding output data. 


Print Position 1 Print Position 67 
eee 1 BLOCK 1 FIELD 1 GROUP 6 BLOCK 6 FIELD é 
7 : (output data) x z a 
GROUP 2 BLOCK 2 FIELD 2 GROUP 7 BLOCK 7 FIELD 7 
: 7 eucoue data) : : 
GROUP 5 BLOCK 5 FIELD 5 GROUP 0 BLOCK O FIELD 0 


— b aaad = —_ —_ — 


= = (output data) = = 7 


Literal Data In a Format Statement 


Literal data consists of a string of alphameric and special charac- 
ters written within the FORMAT statement and enclosed in apostrophes. 
The string of characters must be less than or equal to 255. For 
example: 


25 FORMAT (' 1964 INVENTORY REPORT") 


An apostrophe character within the string is represented by two 
successive apostrophes (either with or without embedded blanks). For 
example, the characters DON'T are represented as: 


DON‘ *T 


The effect of the literal format code depends on whether it is used 
with an input or output statement. 


INPUT 


A number of characters, equal to the number of characters between the 
apostrophes, are read from the designated data set. These characters 
replace, in storage, the characters within the apostrophes. 


For example, the statements: 


5 FORMAT (° HEADINGS‘) 


READ (3,5) 


would cause the next 9 characters to be read from the data set 
associated with data set reference number 3; these characters would 
replace the blank and the 8 characters H,E,A,D,i,N,G, and S in storage. 
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OUTPUT 


All characters (including blanks) within the apostrophes are written 
as part of the output data. Thus, the statements: 


5 FORMAT (° THIS IS ALPHAMERIC DATA") 


WRITE (2,5) 


would cause the following record to be written on the data set 
associated with the data set reference number 2: 


THIS IS ALPHAMERIC DATA 


H Format Code 


a a a eH 


| General Form 


Where: wis an unsigned integer constant that is less than or equal 


to 255 and specifies the number of characters following H. 


The H format code is’ used in conjunction with the transferral of 
literal data. 


The format code wH is followed in the FORMAT statement by w 
characters. For example, | 


5 FORMAT (31H THIS IS ALPHAMERIC INFORMATION) 


Blanks are significant and must be included as part of the count w. 
The effect of wH depends on whether it is used with input or output. 


1. On input, w characters are extracted from the input record and 
replace the w characters of the literal data in the FORMAT 
statement. 

2. On output, the w characters following the format code are written 
as part of the output record. 
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X Format Code 


LD AEA CED SSA DES ATER PEE MD EOE REREAD GD AED TOLD CRD OAS ARES OES ERAS ERED COE MONE WE COED SEER ENED EES COS EEE OEE ES SE CEE SORE GRUB ATE REED CEES GERD CHET EPRD CREE GREED GAY OED STS GED AD TED EURO CORT ET EEE EE GE GD GE! SE ED ED EES Cie aE ee oe oe I 


| General Form 


wx 


to 255 and specifies the number of blanks to be inserted on 
output or the number of characters to be skipped on input. 


| 
| 
| Where: wis an unsigned integer constant that is less than or equal 
| 
| 
L 


fee cman: cues quem eee somes qn aul coe 


When the wxX format code is used with a READ statement (i.e., on 
input), w characters are skipped before the data is read in. For 
example, if a card has six 10-column fields of integer quantities, and 
it is not desired to read the second quantity, then the statement: 


5 FORMAT (110,10X,4110) 
may be used, along with the appropriate READ statement. 


When the wxX format code is used with a WRITE statement (i.e., on 
output), w characters are filled with blanks. Thus, the facility for 
Spacing within a printed line is available. For example, the statement: 


(10 FORMAT ('x', 3(F6.2,5X)) 


may be used with an appropriate WRITE statement to print a line as 
follows: 


123. 45bbbbb817. 32bbbbb524.67bbbkb 


T Format Code 


AP ED TD AEE TD ERD AONE EER EE EE AD AD AND A ED AD OD A AD A A A AS TS A AS AD NS ANS NS NE OE AS NN TNS ES ERD ED ON SY EE NY SES TEED GES ED CYS NEED GED NERY GENET RENN CRED ED SND GRIND SED GENO GENE GRO GENID URED COUNTS GHEE URED Gu <ueE cme woe 


to 255 and specifies the position in a FORTRAN record where 


| 
: 
| Where: wis an unsigned integer constant that is less than or equal 
| 
| the transfer of data is to begin. 


Input and output may begin at any position by using the format code 
Tw. Only when the output is printed does the correspondence between w 
and the actual print position differ. In this case, because of the 
carriage control character, the print position corresponds to w-1, as 
may be seen in the following example: 


5 FORMAT (T40, £1964 INVENTORY REPORT', T80, ‘DECEMBER‘', Ti, ‘° PART 
1 NO. 10095") 
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The preceding FORMAT statement would result in a printed line as 
follows: 


Print Print Print 
Position 1 Position 39 Position 79 
t sf 4 

PART NO. 10095 1964 INVENTORY REPORT DECEMBER 


The following statements: 


5 FORMAT (T40, * HEADINGS") 


READ (5,5) 


would cause the first 39 characters of the input data to be skipped, and 
the next 9 characters would then replace the blank and the characters 
H,E,A,D,1,N,G and S in storage. 


The T format code may be used in a FORMAT statement with any type of 
format code. For example, the following statement is valid: 


5 FORMAT (T100, F10.3, T50, E9.3, T1, * ANSWER IS‘) 


Scale Factor - P 


The representation of the data, internally or externally, may be 
modified by the use of a scale factor followed by the letter P preceding 
a format code. The scale factor must be less than or equal to 127 and 
is defined as follows: 


; ; scale factor 
external quantity = internal quantity x 10 


~ INPUT 

For input, when scale factors are used in a FORMAT statement, they 
have effect only on real or double precision data which does not contain 
an E or D format code, respectively. For example, if input data is in 
the form xx.xxxx and it is desired to use it internally in the form 
«XXXxxx, then the format code used to effect this change is 2PF7.4. 

As another example, consider the following input data: 

27 bbb-93.2094bb-175. 8041 bbbb55. 3647 

where b represents a blank. 


The following statements: 


5 FORMAT (12,3F11.4) 


READ (6,5) K,A,B,C 
would cause the variables in the list to assume the following values: 
K : 27 B : ~175.8041 
A 3: -93.2094 C : 55.3647 
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The following statements: 


5 FORMAT (12,1P3F11.4) 


READ (6,5) K,A,B,C 
would cause the variables in the list to assume the following values: 


K : 27 B 
A: -9.3209 C 


~17.5804 
5.5364 


The following statements: 


5 FORMAT (12,-1P3F11.4) 


READ (6,5) K,A,B,C 
would cause the variable in the list to assume the following values: 


K : 27 B : ~1758.041x 
A: -932.094x C : 553.647x 


where the x represents an extraneous digit. 


OUTPUT 
Assume that the variables K,A,B, and C have the following values: 


K : 27 B : -175.8041 
A : -93.2094 C 3: 55.3647 


then the following statements: 


5 FORMAT (12,1P3F11.4) 


WRITE (4,5) K,A,B,C 
would cause the variables in the list to output the following values: 


K 3: 27 B : ~1758.041x 
A: —-932.094x C : 553.647x 


where the x represents an extraneous digit. 


The following statements: 


5 FORMAT (12,-1P3F11.4) 


WRITE (4,5) K,A,B,C 
would cause the variables in the list to output the following values: 
Ke i 27 B : ~17.5804 
A : -9.3209 C 5.5365 
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For output, when scale factors are used, they have effect only on real 
or double precision data. However, this real or double precision data 
May contain an E or D decimal exponent, respectively. A positive scale 
factor used with real or double precision data which contains an E or D 
decimal exponent increases the number and decreases the exponent. Thus, 
if the real data were in a form using an E decimal exponent, and the 
statement FORMAT (1X,1I2,3E13.3) used with an appropriate WRITE statement 
resulted in the following printed line: 


27bbb-0.932Eb02 bbb-0.175Eb03 bbbb0.553Eb02 


Then the statement FORMAT (1X,1I2,1P3E13.3) used with the same WRITE 
statement results in the following printed output: 


27bbb-9. 320EbD01bbb-1. 758Eb02 bbbb5.536EbD01 


The statement FORMAT (1X,12,-1P3E13.3) used with the same WRITE state- 
ment results in the following printed output: 


27bbb-0 .093Eb03bbb-0.017Eb04 bbbbO.055Eb03 


The scale factor is assumed to be zero if no other valiue has been 
given. However, once a value has been given, it will hold for all 
format codes following the scale factor within the same FORMAT state- 
ment. This also applies to format codes enclosed within an additional 
pair of parentheses. Once the scale factor has been given, a subsequent 
scale factor of zero in the same FORMAT statement must be specified by 
OP. 


Carriage Control 


When records written under format control are prepared for printing, 
the following convention for carriage control applies: 


First Character Carriage Advance Before Printing 


Blank One Line 

0 Two lines 

1 To first line of the next page 
+ No advance 


The first character of the output record may be used for carriage 
control and is not printed. It appears in all other media as data. 


Carriage control can be specified in either of two forms of literal 
data. The following statements are equivalent and each would cause two 
lines to be skipped before printing: 

10 FORMAT (0°, 5(F7.3)) 


10 FORMAT (1HO, 5(F7.3)) 
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END FILE STATEMENT 


CP ag ee pe ek ee eS eg ee eg eee ee SE Eg, ee Ce a Se ee ee 


| General Form ] 


}---------------~---~--------------------------------------------------] 


Where: ais an unsigned integer constant or integer variable that 


represents a data set reference number. 


The END FILE statement defines the end of the data set associated 
with a. A subsequent WRITE statement defines the beginning of a new 
data set. 


REWIND STATEMENT 


Where: a is an unsigned integer constant or integer variable that 
represents a data set reference number. 


The REWIND statement causes a subsequent READ or WRITE statement 
referring to a to read data from or write data into the first data set 
associated with a. 


BACKSPACE STATEMENT 


ES RET NED CRE ND COUNED GEER QUENT ED TTD CERN ERD EE CEU RD ATED NTIS AUTRE EEE ACUED OES AEE GOED EETED CREED GEUNE GRAND ERED CURED GED GED ETD ERED EY ED EEE EE EEE TED GERD RED END AEE ENED QUT GNU TEED GREED ERAS GND CUNES GEED GUND EY EAD OEE SEND GERD GENE GED GREE GUND CED SEED QUIRED GEE sta eee 


BACKSPACE a 


Where: a is an unsigned integer constant or integer variable that 
represents a data set reference number. 


The BACKSPACE statement causes the data set associated with a to 
backspace one record. If the data set associated with a is already at 
its beginning, execution of this statement has no effect. 
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DIRECT ACCESS INPUT/OUTPUT STATEMENTS 


There are four direct access I/O statements:1 FILE, and FIND. The 
READ and WRITE statements cause transfer of data into or from internal 
storage. These statements allow the user to specify the location within 
a data set from which data is to be read or into which data is to be 
written. 


The DEFINE FILE statement specifies the characteristics of the data 
set(s) to be used during a direct access operation. The FIND statement 
overlaps record retrieval from a direct access device with computation 
in the program. In addition to these four statements, the FORMAT 
Statement (described previously) specifies the form in which data is to 
be transmitted. The direct access READ and WRITE statements, and the 
FIND statement are the only I/O statements that may refer to a data set 
reference number defined by a DEFINE FILE statement. 


DEFINE FILE Statement 


The DEFINE FILE statement is a specification statement that describes 
the characteristics of any data set to be used during a direct access 
input/output operation. To use the direct access READ, WRITE, and FIND 
statements in a program, the data set(s) must be described with a DEFINE 
FILE statement. Each data set must be described once, and this 
description may appear once in each program or subprogranm. 


Because the DEFINE FILE statement is a specification statement, it 
must precede the first executable statement of the source program as 
well as all statement function definition statements. The description 
must appear logically before the use of an input/output statement with 
the same data set reference number; subsequent descriptions have no 
effect. 


Example: 
DEFINE FILE 2(50,100,L, 12), 3(100,50,L,J3) 


This DEFINE FILE statement describes two data sets, referred to by 
data set reference numbers 2 and 3. The data in the first data set 
consists of 50 records, each with a maximum length of 100 storage 
locations. The L specifies that the data is to be transmitted either 
with or without format control. I2 is the associated variable that 
serves aS a pointer to the next record. 


The data in the second data set consists of 100 records, each with a 
Maximum length of 50 storage locations. The L specifies that the data 
is to be transmitted either with or without format control. J3 is the 
associated variable that serves as a pointer to the next record. 


If an E is substituted for the L in the preceding DEFINE FILE 
Statement, a FORMAT statement is required and the data is transmitted 
under format control. If the data is to be transmitted without format 
control, the DEFINE FILE statement can be written as: 


DEFINE FILE 2(50, 25,0, 12),3(100,13,U,J73) 


EPP ARP RD CATED CIEE AED RECEDED AEE TED END GENET AEE RD END NED ATED SEP OED GD 


1The direct access I/O statements are not available in Basic Programming 
Support FORTRAN IV. 
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| General Form 


DEFINE FILE a, (m4 ,¥s,f1,V1) «a2 (Ma,Fa,fasVa),~---.A8n(Mn, Ins Ens Vn) 


Where: a represents an unsigned integer constant that is the data 
set reference number. 


m represents an unsigned integer constant that specifies the 
number of records in the data set associated with a. 


Xr represents an unsigned integer constant that specifies the 
maximum size of each record associated with a. The record 
size is measured in characters, storage locations, or 
storage units. (A storage unit is the number of storage 
locations divided by four and rounded to the next highest 
integer). The method used to measure the record size 
depends upon the specification for f. 


£ specifies that the data set is to be read or written 
either with or without format control; £ may be one of the 
following letters: 


L indicates that the data set is to be read or written 
either with or without format control. The maximum record 
size is measured in number of storage locations. 


E indicates that the data set is to be read or written 
under format control (as specified by a FORMAT statement). 
The maximum record size is measured in number of charac- 
ters. 


U indicates that the data set is to be read or written 
without format control. The maximum record size is 
measured in number of storage units. 


v represents a nonsubscripted integer variable called an 
associated variable. At the conclusion of each read or 
write operation, v is set to a value that points to _ the 
record that immediately follows the last record transmitted. 
At the conclusion of a find operation, v is set to a value 
that points to the record found. 


Re SARS ER ear eR AID em WOES CURRED CGR «RE SERED CRATE GCOS CTE EGE MND ARE QRS ED GUE ORIYA CEGED SETEED CEPER REDD “OIE comterD <CERED EMIT SREOR GREER SESE CREED CRED GESTED GRUNGE CORED SEED “GRINS NWA GETEIED CRD RE EEE COREE GENTS TA ERED SOE REED GENES GREED GED RP GET SEEETS GOED cAtHR “OEE CURED GER ED SE ci GEES aEE GO eR 


f° mee eee eee cme cee ee ce ee ee ee ee ee ee oe oe coe ee ee ee ee ee ee ee ee eee ee ED ecm seen Se ee ee comes eee eee ee cee ee eee eee ey 


Programming Considerations 


When programming for direct access input/output operations, the user 
must establish a correspondence between FORTRAN records and the records 
described by the DEFINE FILE statement. All of the conventions of 
FORMAT control discussed in the section “FORMAT STATEMENT" are applica- 
ble. 


For example, to process the data set described by the statement: 
DEFINE FILE 8(10,48,L,K8) 


the FORMAT statement used to control the reading or writing could not 
specify a record longer than 48 characters. The statements: 


FORMAT (4F12.1) or 
FORMAT (112,9F4.2) 
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define a FORTRAN record that corresponds to those records described by 
the DEFINE FILE statement. The records can also be transmitted under 
FORMAT control by substituting an E for the L and rewriting the DEFINE 
FILE statement as: 


DEFINE FILE 8(10,48,E,K8) 


Programs may share an associated variable only as a COMMON variable. 
The following example shows how this can be accomplished. 


COMMON. IUAR SUBROUTINE SUBI(A,B) 
DEFINE FILE 3(100,10,L, IUAR) COMMON IUAR 


ITEMP=IUAR 

CALL SUBI(ANS, ARG) 
8 IF (IUAR-ITEMP) 20,16, 20 
20 CONTINUE 


In this example, the program and the subprogram share the associated 
variable IUAR. An input/output operation that refers to data associated 
with data set reference number 3 and is performed in the subroutine 
causes the value of the associated variable to be changed. The 
associated variable is then tested in the main program in statement 8. 


READ Statement 


The READ statement causes data to be transferred from a direct access 
device into internal storage. The data set being read must be defined 
with a DEFINE FILE statement. 


| General Form 


AEs IE NEE ENS EEN SINED IY EO AS A ES ID ACN EI ON OED IS AS ES TE SEED ENE ND EE NY ILE SND CANE GER GURNEE CN AE ERED SUED CURRED ATE RN CORED GUND GD AEN CGT RED SNE EAR EEIT SS GRINS ETRE SEY GNU ESE REINS SERIE SERED ENE AEE IETS MIRREN ET GNM REND CEO GEE QD OREN meT <Ph EAD 


READ (a‘xr,b) list 


Where: a iS an unsigned integer constant or integer variable that 
represents a data set reference number; a must be followed 
by an apostrophe ('). 


rk is an integer expression that represents the relative 
position of a record within the data set associated with a. 


b is optional and, if given, is the statement number of the 
FORMAT statement that describes the data being read. 


list is a series of variable or array names, which may be 
indexed and must be separated by commas. These names 
specify the number of items to be read and the storage 
locations into which the data is to be placed. The list has 
the same forms and conventions as the list for the sequen- 
tial READ statements. 


ND Ge CELE SANT FNS EGND CREP CERES COTE EAE SLED CNT CES ED CORIO ETS GP AIUD ALTE ANS TID OME IUD ETS CEILI RING AED SNIED CUED EEE GG GED GUND GERD CRED ESD IED GTS GU CTD TINS CENUD EdCRD GENTE CEES CUD OND ARTE GHED KOEED ERED CON ANID GEARY THUD GANS QC GOED ND CHINN TINY CED NONE SID CMS CON CR eee OTE 


FFP SEED ee Stee See Cee SEEN Gee SED ER Clee ee i Se a ee cee ates Sp 
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Example: 


DEFINE FILE 1(500,100,L,1ID1) ,2(100, 28,L,ID2) 
DIMENSION M(10) 


10 FORMAT (5120) 
9 READ (1'16,10) (M(K),K=1,10) 


READ statement 9 transmits data from the data set associated with 
data set reference number 1, under control of FORMAT statement 10; 
transmission begins with record 16. Ten data items of 20 characters 
each are read as specified by the I/O list and FORMAT statement 10. Two 
records are read to satisfy the I/O list, because each record contains 
only five data items (100 characters). The associated variable ID1 is 
set to a value of 18 at the conclusion of the operation. 


READ statement 13 transmits data from the data set associated with 
data set reference number 2, without format control; transmission begins 
with record 26. Data is read until the I/O list for statement 13 is 
satisfied. Because the DEFINE FILE statement for data set 2 specified 
the record length as 28 storage locations, the I/O list of statement 13 
calls for the same amount of data (the seven variables are type real and 
each occupies four storage locations). The associated variable ID2 is 
set to a value of 27 at the conclusion of the operation. If the value 
of ID2 is unchanged, the next execution of statement 13 reads record 32. 


The DEFINE FILE statement in the previous example can also be written 
as 
DEFINE FILE 1(500,100,E,1ID1),2(100,7,U,ID2) 
The FORMAT statement may also control the point at which reading 
starts. For example, if statement 10 in the example was 
10 FORMAT (7/5120) 
records 16 and 17 are skipped, records 18 and 19 are read, and ID1 is 


set to a value of 20 at the conclusion of the read operation in 
Statement 9. 


WRITE Statement 


The WRITE statement causes data to be transferred from internal 
storage to a direct access device. The data set being written must be 
defined with a DEFINE FILE statement. 
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WRITE (atry,b) list 


Where: ais an unsigned integer constant or integer variable that 
represents a data set reference number; a must be followed 
by an apostrophe ('‘'). 


r is an integer expression that represents the relative 
position of a record within the data set associated with a. 


FORMAT statement that describes the data being written. 


list is a series of variable or array names, which may be 
indexed and must be separated by commas. These names 
Specify the number of items to be written and the locations 
in storage from which the tata is to he taken. The iist has 
the same forms and conventions as the I/O list for the 


sequential WRITE statements. 


URED QUES EES ED ENT ERE GE ES STEED ACKER ATR OEE ETRE SOS ATE GEER CEES ATTEN GENE ADIGE COUKRY SURED etunie SEED SOUND GREETS GREER RIED CIEIOU AERA EES, 02111) USE MLCT CIE ORL ANN TTD AID IONS TSO RECT ET RRR COTE CUTER CITES CLS ROOE SIO ABNEY ZIG CHTTRD RIERA AIST EI (AONE ACRE SAIRRID ED ORO OE EEDLTD ARENDS SAID ED ape SEI 
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Example: 


DEFINE FILE 1(500,190,L, 201), 2¢(106, 28, L, ID2) 
DIMENSION M(10) 
ID2=21 


10 FORMAT (5120) 
8 WRITE (1°16,10) (M(K),K=1,10) 


11 WRITE (2"ID2+5) A,B,C,D,E,F,G 


WRITE statement 8 transmits data into the data set associated with 
the data set reference number 1, under control of FORMAT statement 10; 
transmission begins with record 16. Ten data items of 20 characters 
each are written as specified by the I/O list and FORMAT statement 10. 
Two records are written to satisfy the I/O list because each record 
contains 5 data items (100 characters). The associated variable ID1 is 
set to a value of 18 at the conclusion of the operation. 


WRITE statement 11 transmits data into the data set associated with 
data set reference number 2, without format control; transmission begins 
with record 26. The contents of 28 storage locations are written as 
specified by the I/O list for statement 11. The associated variable ID2 
is set to a value of 27 at the conclusion of the operation. Note the 
correspondence between the records described (28 storage locations per 
record) and the number of items called for by the I/O list (7 variables, 
type real, each occupying four storage locations). 


The DEFINE FILE statement in the previous example can also be written 
as 

DEFINE FILE 1(500,100,E,ID1) ,2(100,7,U,ID2) 

As with the READ statement, a FORMAT statement may also be used to 
control the point at which writing begins. 
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FIND Statement 


The FIND statement permits record retrieval to proceed concurrently 
with computation. By using the FIND statement, the user can increase 
the object program execution speed. There is no advantage to ae FIND 
statement preceding a WRITE statement. The data set from which a record 
is being retrieved must be defined with a DEFINE FILE statement. 


Ee er nc ce ep a a er er nd Fe ae oe na ge ee eo OT CE etre gems gO mE Oe Oe ee 
| General Form | 
}----------------------------------------------------------------- { 
| | 
| FIND (a*r) | 
| | 
| Where: a is an unsigned integer constant or integer variable that | 
| represents a data set reference number; a must be followed | 
i by an apostrophe (‘). | 
| | 
| x is an integer expression that represents the relative | 
i position of a record within the data set associated with a. | 
Ue at ka ere J 
Example: 


10 FIND (3°50) 


15 READ (3°50) A,B 


While the statements between statements 10 and 15 are executed, 
record 50, in the data set associated with data set reference number 3, 
is retrieved. If a WRITE statement refers to this record between the 
issuing of the FIND statement and the READ statement, the FIND operation 
is nullified. 


General Examples -- Direct Access Operations 


Example 1: 


DEFINE FILE 8(1000,72,L,ID8) 
DIMENSION A(100),B(100) ,C (100) ,D(100),E(100),F(100) 


15 FORMAT (6F12.4) 
FIND (8'5) 


ID8=1 

DO 100 I=1,100 

READ (8°ID8+4,15)A(I),B(1I),C(I),D(I) ,E(I), F(T) 
100 CONTINUE 


DO 200 I=1,100 
WRITE (8*°ID8+4,15)A(I), B(I),C(I) ,D(I) ,E(I), F(T) 
200 CONTINUE 


END 
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Explanation: 


Example 1 illustrates the ability of direct access statements to 
gather and disperse data in an order designated by the user. The first 
DO loop in the example fills arrays A through F with data from the 
fifth, tenth, fifteenth,..., and five-hundredth record associated with 
data set reference number 8. Array A receives the first value in every 
fifth record, B the second value and so on, as specified by FORMAT 
Statement 15 and the I/O list of the READ statement. At the end of the 
READ operation, each record has been dispersed into arrays A through F. 
At the conclusion of the first DO loop, ID8 has a value of 501. 


The second DO loop in the example groups the data items from each 
array, aS specified by the I/O list of the WRITE statement and FORMAT 
Statement 15. Each group of data items is placed in the data set 
associated with data set reference number 8. Writing begins at the 
505th record and continues at intervals of five, until record 1000 is 
written. 


Example 2: 
C MAIN PROGRAM C SUBROUTINE ONE 
COMMON I SUBROUTINE SUB1 (AA) 
, COMMON J 
‘ DEFINE FILE 9(100,100,E,9) 
1 READ (1,2) I ‘ 
2 FORMAT (14) ‘ 
I=IABS(T) RETURN 
IF (I) 10,20,10 END 
10 CALL SUB1 (A) 
GO TO 70 


20 CALL SUB2 (A) 
70 CONTINUE 


; C SUBROUTINE TWO 
; SUBROUTINE SUB2 (BB) 
; COMMON K 
DEFINE FILE 9(125,80,L,K) 
WRITE (9°I+1,100) X,Y,Z : 
100 FORMAT (3F10. 3) - 
, RETURN 
: END 
END 


Explanation: 


Example 2 illustrates the use of two different DEFINE FILE statements 
to describe the characteristics of the data set associated with data set 
reference number 9. If subroutine SUB1 is called, the data set contains 
100 records, each with a maximum length of 100 characters; the data is 
to be transmitted under format control, and the associated variable is 
‘a If subroutine SUB2 is called, the data set contains 125 records, 
each with a maximum length of 80 storage locations; the data is to be 
transmitted either with or without format control, and the associated 
variable is K. Because the associated variables are declared to be in 
COMMON along with I, the information is shared between the main program 
and the two subroutines. 
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SPECIFICATION STATEMENTS 


The specification statements provide the compiler with information 
about the nature of the data used in the source program. In addition, 
they supply the information required to allocate locations in storage 
for this data. Specification statements describing data must precede 
any statements which refer to that data. For example, if an element of 
an array is to be made equivalent to a variable, the specification 
Statement that declares the size of the array (e.g., a DIMENSION 
statement) must precede the EQUIVALENCE statement. 


All specification statements must precede the first executable 
statement of the source program.+ They must also precede all Statement 
Function definition statements. Therefore, the source program layout is 
as follows: 


1. Specification Statements. 
2. Statement Function Definition Statements. 
3. Executable Statements. 


FORMAT statements may appear anywhere in the program. They do not 
affect the sequence of execution. 


Explicit Specification Statements 


PG AON CC LER NEY IT AEN tio CAD COED CEE CEP OSG GA ACRES AEE EE A A SED OED OE GID GD ED SO GE TD EES AEE eee CED REED GERNOT GU RD TD ACRE ED ORD CRED FO GE GEES RE GED ECO EE EP GET EEE ACME GSR GESEP GED OEE AUEE ae aiieD GE ED a a GED ee Ee ce 


7 
| General Form 


Type a(k,),b(k2), 2@ o,z(Kn) 
Where: Type is INTEGER, REAJ.. or DOUBLE PRECISION. 


a,b,..e,zZ represent variable, array, or function names (see 
the section, “SUBPROGRAMS") 


(ki), (ka),---, (Kn) are optional. Each k is composed of 1 
through 3 unsigned integer constants, separated by commas, 
representing the maximum value of each subscript in the 
array. 


NS ED ES ED EER TS SOD RC) ED ARENAS EET SE ITE RONNIE ERED GENYT SANE CUED IEE <A CREED TEES CONES EARS SEEMED REDD SEE GILES ORCL ERED CORD GERIED GOURD EUED GUTAD CRAP GUC OGREIRD YEE SEUNEND “RIND SEW AERO CIS NYS ND <I GRUNT RN GUND SERN EST GUD CARED GENE CLEP QUYD QUEERS CORED GUERD CN GERD CRED GED GEE RED auEED aE Ee aoe 


GF see eee we oe ee ee ee ere ee ee eee ee ep 


The Explicit specification statements declare the type (INTEGER, 
REAL, or DOUBLE PRECISION) of a particular variable or array by its 
name, rather than by its initial character. This differs from the other 
way of specifying the type of a variable or array (1.e., predefined 
convention). In addition, the information necessary to allocate storage 
for arrays (dimension information) may be included within the statement. 
However, if this information does not appear in an Explicit specifi- 
cation statement, it must appear in a DIMENSION or COMMON statement 
(see, “DIMENSION Statement" or “COMMON Statement"). 


4In Operating System FORTRAN IV (E) and Basic Programming Support 
FORTRAN IV, EQUIVALENCE statements must follow any DIMENSION, COMMON, or 
explicit specification statements; they need not follow DEFINE FILE 
statements in FORTRAN IV (E). 
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Example 1: 
INTEGER DEV, ARRAY, SMALL 


Explanation: 


This statement declares the type of the variables DEV, ARRAY, and 
SMALL as integer and thus overrides the implied declaration made by the 
predefined convention. 


Example 2: 


REAL ITA,JOB, MATRIX (5,2,6) 


Explanation: 


This statement declares the type of the array, MATRIX, and variables, 
ITA and JOB, as real. In addition, it declares the size (dimension) of 
the array MATRIX. This statement overrides the implied declaration made 
by the predefined convention. 


Example 3: 


DOUBLE PRECISION DOUB, TWIN 


Explanation: 


This statement declares the type of the variables, DOUB and TWIN , as 
double precision. 


DIMENSION Statement 


| General Form 


DIMENSION ay (ky) ,a2 (ka), a3 (K3) 4-20 ¢An (Kn) 
Where: ayy A@ay Azzeeey An AYE array names. 


Kae Kae Kay,----,Kn are each composed of 1 through 3 unsigned 
integer constants, separated by commas, representing the 
maximum value of each subscript in the array. 


[Fe nee ee ES ee Ree atm Ae ates Ee 


The information necessary to allocate storage for arrays used in the 
source program may be provided by the DIMENSION statement. Allocation 
information should be given to an array on its first appearance in a 
source program; however, for subprograms, the SUBROUTINE or FUNCTION 
statement may include a dummy argument that is dimensioned later. The 
following examples illustrate how this information may be declared. 


Examples: 


DIMENSION A (10), ARRAY (5,5,5), LIST (10,100) 
DIMENSION B(25,50), TABLE(5,8,4) 
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COMMON Statement 


Bee ee eg he a, Ve ete PAGE ge ee pepe Va Fee me Oe TN eR Gee ey ee ape te ee 


| General Form i 


}--~—---- ~~ === 2 nn nnn nnn 


COMMON ag (ki ),a2 (Kka),a3 (ks) peee gan (Ky) 


Where: €@4,82,83,---yay are variable or array names. 


through 3 unsigned integer constants, separated by commas, 
representing the maximum value of each subscript in the 


[ 
| 
| 
i KaeKa,Kay---,Kn are optional and are each composed of 1 
| 
| 
| array. 


| 


The COMMON statement may be used to provide dimension information. 
However, this information must be given on the first appearance of an 
array. For example, if an array has appeared first in an Explicit 
specification statement, the array should be dimensioned in that 
statement rather than in the COMMON statement. 


Variables or arrays that appear ina calling program or a subprogram 
may be made to share the same storage locations with variables or arrays 
in other subprograms by use of the COMMON statement. For example, if 
one program contains the statement: 

COMMON TABLE 
and a second program contains the statement: 
COMMON TREE 
the variable names TABLE and TREE refer to the same storage locations. 


If the main program contains the statements: 


REAL A,B,C 
COMMON A,B,C 


and a subprogram contains the statements: 


REAL X,Y,Z 
COMMON %,Y,2 


then A shares the same storage location as X, B shares the same storage 
location as Y¥, and C shares the same storage location as Z. 


Common entries appearing in COMMON statements are cumulative in the 
given order throughout the program; that is, they are cumulative in the 
sequence in which they appear in all COMMON statements. For example, 
consider the following two COMMON statements: 


COMMON A, B, C 
COMMON G, H 


These two statements have the same effect as the single statement: 


COMMON A, By Cy G, H 
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Redundant entries are not allowed. For example, the following 
Statement is invalid: 


COMMON A,B,C,A 


Consider the following examples: 


Example 1: 
Calling Program Subprogram 
3 SUBROUTINE MAPMY (...) 


COMMON A, B, C, R(100) 


REAL A,B,C COMMON X%, ¥, Z, S(€100) 
INTEGER R REAL X,Y¥,24 

. INTEGER S 
CALL MAPMY (...) « 


Explanation: 


In the calling program, the statement COMMON A,B,C,R(100) would cause 
412 storage locations (4 locations per variable) to be reserved in the 
following order: 


Beginning of COMMON area ,—-—--~—--—----—--- 7 
| A | 
[--------------- j 
B | 
}--------------- 4 
| Cc | 
|--------------- { 
i R(1) i 
| R(2) i 
. | 
° | 
| ° | 
| R(100) | 
ee eee ee ae eee 


4 storage locations 


The statement COMMON X, Y, 4, S(100) would then cause the variables 
X, Y, Z, and S(1i)...S(100) to share the same storage space as A, B, C, 
and R(1)...R(100), respectively. 


From the above example, it can be seen that COMMON statements may be 
used to serve an important function: namely, as a medium to implicitly 
transmit data from the calling program to the subprogram. That is, 
values for X, Y, Z, and S(1)...S(100), because they occupy the same 
storage locations as A, B, C, and R(i)...R(100), do not have to he 
transmitted in the argument list of a CALL statement. Arguments passed 
through COMMON must follow the same rules of presentation with regard to 
type, length, etc., as arguments passed in alist. (See the section, 
"“SUBPROGRAMS.") 
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Example 2: 


Assume COMMON is defined ina main program and 3 subprograms’ as 
follows: 


Main program: COMMON A,B,C, (A and B are 8 storage locations, 
C is 4 storage locations) 
Subprogram 1: COMMON D,E,F (D and E are 8 storage locations, 
F is 4 storage locations) 
Subprogram 2: COMMON Q,R,S,T,U (4 storage locations each) 
Subprogram 3: COMMON V,W,X,Y,2 (4 storage locations each) 


The correspondence of these variables within COMMON can be illustrated 
as follows: 


Main Program Subprogram 1 Subprogram 2 Subprogram 3 
COMMON A,B,C COMMON D,E,F COMMON Q,R,S,T,U COMMON V,W,X,Y,2Z 
ore 1 Pernt eo 1 ES eer 1 laeia beara rG aati: 1 
| | | | Q paar | V | 
[a2 Shs eS De }---------- { b--------—- 
| | | | R +» W | 
}-~-------- { |---------- { |~--------- { }---------- 
| | | | S Ss X | 
[Se Sees ee | }-----~---- { aaa 
| | | | T ~—____—_—>| ¥ | 
[------~-~- { [-----—-~--4 [---------- { }---—------ { 
| C fge-_-__—-+ F fg Uk +} Z| 
Ue ceeete es SS, om Wieder nea eee ee J Ue ce J Seer een er J 
4 storage 4 storage 4& storage 4 storage 
locations locations | locations locations 


In this case, the variables A,B,C and D,E,F may be validly referred to 
in their respective programs as may Q,R,S,T,U and V,W,X,Y,Z. In 
addition, Subprogram 1 may implicitly refer to C,U, and Z by explicitly 
referring to F. 


To insure proper boundary alignment, the user is advised to arrange 
the variables in COMMON in the following order: 


Double Precision 
Real or Integer 


EQUIVALENCE Statement 


| General Form 


EQUIVALENCE (a, b, c, «--), (d, e, f£,--.) 


Where: a, b, c, qd, e, £,--- are variables that may be subscripted. 
The subscripts may have two forms: If the variable is singly 
subscripted it refers to the position of the variable in the 
array (i.e., first variable, 25th variable, etc). If the 
variable is muiti-subscripted, it refers to the position in 
the array in the same fashion as the position is referred to 
in an arithmetic statement. 


(awe om coe Oe ee ee ee ee oe ee 


The EQUIVALENCE statement provides the option for controlling the 
allocation of data storage within a single program or subprogram. It is 
analogous to the option of using the COMMON statement to control the 
allocation of data storage among several programs. In particular, when 
the logic of the program permits it, the number of storage locations 
used can be reduced by causing locations to be shared by two or more 
variables of the same or different types. 


In storage established by EQUIVALENCE statements, double precision 
variables must precede real or integer variables. 


Example 1: 


DIMENSION B(5), C(10, 10), D(5, 10, 15) 
EQUIVALENCE (A, B(1), C(5,3)), (D(5,10,2), E) 


Explanation: 


This EQUIVALENCE statement indicates that the variables A, B(1), and 
C(5,3) are to share the same storage locations. All other corresponding 
variables in arrays B and C share the same storage locations; e.g., B(2) 
and C(6,3). In addition, it specifies that D(5,10,2) and E are to share 
the same storage locations. In this case, the subscripted variables 
refer to the position in an array in the same fashion as the position is 
referred to in an arithmetic statement. Note that variables or arrays 
that are not mentioned in an EQUIVALENCE statement are assigned unigue 
Storage locations. The EQUIVALENCE statement must not contradict itself 
Or any previously established equivalences. For example, the further 
equivalence specification of B(3) with any other element of the array C, 
other than C(7,3), is invalid. 


Example 2: 


DIMENSION B(5), C(10, 10), D(5, 10, 15) 
EQUIVALENCE (A, B(1), C(25)), (D(100), E) 


Explanation: 


This EQUIVALENCE statement indicates that the variable A, the first 
variable in the array B, namely B(1), and the 25th variable in the array 
C, namely C(5,3), are to share the same storage locations. In addition, 
it also specifies that D(100) (i.e., D(5,10,2) ) and E are to share the 
Same storage locations. Note that the effect of the EQUIVALENCE 
statement in Examples 1 and 2 is the same. 


Variables that are brought into COMMON through EQUIVALENCE statements 
May increase the size of the block as indicated by the following: 
Statements: 


COMMON A, B, C 
DIMENSION D(3) 
EQUIVALENCE (B,D(1)) 


This would cause a common area to be established containing the 
variables A, B, and C. The EQUIVALENCE statement would then cause the 
variable D(1) to share the same storage location as B, D(2) to share the 
same storage location as C, and D(3) would extend the size of the common 
area, in the following manner: 


A (lowest location of the common area) 
B, D(1) 
C, D2) 

D(3) (highest location of the common area) 
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Since arrays are stored in consecutive forward locations, a variable 
may not be made equivalent to another variable of an array in such a way 
as to cause the array to extend before the beginning of the common area. 
For example, the following EQUIVALENCE statement is invalid: 


COMMON A, B, C 
DIMENSION D(3) 
EQUIVALENCE (B, D(3)) 


because it would force D(1) to precede A, as follows: 


D(1) 
A, D(2) (lowest location of the common area) 
B, D(3) 
C (highest location of the common area) 


Two variables in COMMON may not be made equivalent. In addition, a 
user must observe the suggestion on boundary alignment mentioned in the 
description of the COMMON statement. Any double precision variable that 
is made equivalent to any variable in COMMON must be aligned properly. 
For example, the following statements wili produce proper boundary 
alignment: 


REAL A,B 

DOUBLE PRECISION DOUB, TWIN 
INTEGER I,J 

COMMON DOUB,A,I,B 
EQUIVALENCE (A, TWIN), (B,J) 
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SUBPROGRAMS 


It is sometimes desirable to write a program which, at various 
points, requires the same computation to be performed with different 
data for each calculation. It would simplify the writing of that 
program if the statements required to perform the desired computation 
could be written only once and then could be referred to freely, with 
each subsequent reference having the same effect as though these 
instructions were written at the point in the program where the 
reference was made. 


For example, to take the cube root of a number, a program must be 
written with this object in mind. If a general program were written to 
take the cube root of any number, it would be desirable to be able to 
combine that program (or subprogram) with other programs where cube root 
calculations are required. 


The FORTRAN language provides for the above situation through the use 
of subprograms. There are three classes of subprograms: Statement 
Functions, FUNCTION subprograms, and SUBROUTINE subprograms. In addi- 
tion, there is a group of FORTRAN supplied subprograms (see Appendix C). 


The first two classes of subprograms are called functions. Functions 
differ from SUBROUTINE subprograms in that functions return at least one 
value to the calling program; whereas, SUBROUTINE subprograms need not 
return any. 


NAMING SUBPROGRAMS 


A subprogram name consists of from 1 through 6 alphameric characters, 
the first of which must be alphabetic. A subprogram name may not be a 
variable name and may not contain special characters other than the 
blanks (see Appendix A). Blanks embedded in a subprogram name are 
ignored. The type of a subprogram can be indicated in the same manner 
as variables. 


1. Type Declaration of a Statement Function: Such declaration may be 
accomplished in one of two ways: by the predefined convention, or 


by the Explicit specification statements. Thus, the same rules for 
declaring the type of variables apply to Statement Functions. 


2. Type Declaration of FUNCTION Subprograms: The declaration can be 


made in one of two ways: by the predefined convention or by an 
explicit specification (see the section, “Type Specification of the 
FUNCTION Subprogram"). In addition, the type may appear in the 
FUNCTION definition statement. 


3. Type Declaration of a SUBROUTINE Subprogram: The type of a 
SUBROUTINE subprogram can not be defined because the results that 


are returned to the calling program are dependent only on the type 
of the variable names appearing in the argument list of the calling 
program and/or the implicit arguments in COMMON. 
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FUNCTIONS 


A function is a statement of the relationship between a number of 
variables. To use a function in FORTRAN, it is necessary to: 


1. Define the function (i.e., specify what calculations are to be 
performed). 
2. Refer to the function by name where required in the program. 


Function Definition 


There are three steps in the definition of a function in FORTRAN: 


1. The function must be assigned a unique name by which it may be 
called (see the section "Naming Subprograms"). 

2. The arguments of the function must be stated. 

3. The procedure for evaluating the function must be stated. 


Items 2 and 3 are discussed in detail in the sections dealing with 
the specific subprogram (e.g., "Statement Functions", "FUNCTION 
Subprograms", etc.). 


Function Reference 


When the name of a function appears in any FORTRAN arithmetic 
expression, this, effectively, references the function. Thus, the 
appearance of a function with its arguments in parentheses causes the 
computations to be performed as indicated by the function definition. 
The resulting quantity replaces the function reference in the expres- 
Sion. The type of the name used for the reference must agree with the 
type of the name used in the definition. 


STATEMENT FONCTIONS 


Statement functions are defined by a Single arithmetic statement 
within the program in which they appear. For example, the statement: 


FUNC(A,B) = 3. *AtBt¥*¥2.+Xt+Y+Z 


defines the statement function FUNC, where FUNC is the function name and 
A and B are the function arguments. 


The expression on the right defines those computations which are to 
be performed when the function is used in an arithmetic statement. This 
function might be used in a statement as follows: 


C = FUNC(D,E) 
which is equivalent to: 
C = 3.*DtE**¥2.+X+YtZ 
Note the correspondence between A and B in the function definition 
Statement and D and E in the arithmetic statement. The quantities A and 
B enclosed in parentheses following the function name are the arguments 
of the function. They are dummy variables for which the quantities D 


and E, respectively, are punStr eure? when the function is used in an 
arithmetic statement. 
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: 
| General Form [ 


}--------------------------------------------------------------------+4 


name (a,b,...,n) = expression 


Where: name is any subprogram name (see the section “Naming 
Subprograms"™). 


Scripted variables. 


expression iS any arithmetic expression that does not 
contain subscripted variables. Any statement functions 


| 
| | 
| | 
| | 
| | 
| | 
| a,b,...,n are distinct (within the same statement) nonsub- | 
| | 
| | 
| | 
| | 
| appearing in this expression must be defined previously. | 


A maximum of 15 variables appearing in the expression may be used as 
arguments of the function. The actual arguments must correspond in 
order, number, and type to the dummy arguments. There must be at least 
one argument. 


Note: All Statement Function definitions to be used in a program must 
precede the first executable ~*.atement of the program. 


Examples: 


Valid statement function definitions: 


SUM(A,B,C,D) = AtBtCtD 

FUNC(Z) = AtX*Y*Z 

AVG(A,B,C,D) = (AtBt+Ct+D) 74 

ROOT(A, B,C) = SQRT(A**24+B**2+C*#*2) 


Note: The same dummy arguments may be used in more than one Statement 
Function definition and may be used as variables outside Statement 
Function definitions. 





Invalid statement function definitions: 


SUBPRG(3,d, K) =3*1+5**3 (arguments must be variables) 
SOMEF (A(I) , B) =A(I) /Bt+3. (arguments must be nonsub- 
scripted) 


SUBPROGRAM (A, B) =A* *2+B**2 (function name exceeds limit 
of six characters) 


3FUNC (D)=3.14*E (function name must begin with 
an alphabetic character) 

ASF (A) =A+B(T) (subscripted variable in the 
expression) 


GRADE = AVG(ALAB, TERM, SUM(TESTI, TEST2, TEST3, TEST4), FACTOR) 
Invalid statement function references: 


WRONG = SUM(TAX, FICA) (number of arguments 

| does not agree with 
above definition) 

MIX = FUNC(T) (mode of argument 
does not agree with 
above definition) 
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FUNCTION SUBPROGRAMS 


The FUNCTION subprogram is a FORTRAN subprogram consisting of any 
number of statements. It is an independently written program that is 
executed wherever its name appears in another program. 


c 
| General Form 





1 

| 
PT 
| FUNCTION name (ag¢ao 1Asee e °*4an) | 
, | 
= | 
| RETURN | 
| . | 
_—_ | 
| END | 
| Where: name is  subprogram name (see the section "Naming | 
| Subprograms"). | 
| | 
[ Q@4¢42sAgeeee48n are nonsubscripted variable or array names, | 
| or the dummy names of SUBROUTINE or other FUNCTION subpro- | 

rams. (There must be at least one argument in the argument 
Es g 

i list.) | 
a a a a a ee J 


Since the FUNCTION is a separate subprogram, the variables and 
statement numbers within it do not relate to any other program. 


The FUNCTION subprogram may contain any FORTRAN statement except a 
SUBROUTINE statement or another FUNCTION statement. 


The arguments of the FUNCTION subprogram (1.e., a4,a2,83,---,an) may 
be considered to be dummy variable names. These are replaced at the 
time of execution by the actual arguments supplied in the function 
reference in the calling program. The actual arguments may be any of 
the following: 


1. any type of constant 

2. any type of subscripted or nonsubscripted variable 

3. an array name 

4. an arithmetic expression 

5. the name of another FUNCTION or SUBROUTINE subprogram. 


The actual arguments must correspond in number, order, and type to the 
dummy arguments. The array size must also be the same. The name of the 
FUNCTION subprogram cannot be typed with an Explicit specification 
statement in the subprogram. 


The relationship between variable names used as arguments in the 
calling program and the dummy variables used as arguments in the 
FUNCTION subprogram is illustrated in the following example: 
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Example 1: 


Calling Program FUNCTION Subprogram 
‘ FUNCTION SOMEF (X, Y) 
3 SOMEF = X/Y 

A = SOMEF(B,C) RETURN 

= END 


Explanation: 


In the above example, the value of the variable B of the calling 
program is used in the subprogram as the value of the dummy variable xX; 
the value of C is used in place of the dummy variable Y. Thus if B= 
10.0 and C = 5.0, then A = B/C, which is equal to 2.0. 


The name of the function must be assigned a value at least once in 
the subprogram as the argument of a CALL statement, as the variable name 
on the left side of an arithmetic statement, or in an input list (READ 
statement) within the subprogram. 


Example 2: 
Calling Program FUNCTION Subprogram 
FUNCTION CALC (A,B,J) 
é I = J¥2 
ANS = ROOT1*CALC(X, Y, I) : 
: CALC = A*¥*I/B 
RETURN 
END 


Explanation: 


In this example, the values of X, Y, and I are used in the FUNCTION 
Subprogram as the values of A, B, and J, respectively. The value of 
CALC is computed, and this value is returned to the calling program 
where the value of ANS is computed. The variable I in the argument list 
of CALC in the calling program is not the same as the variable I 
appearing in the subprogram. 


When a dummy argument is an array name, an appropriate DIMENSION or 
Explicit specification statement must appear in the FUNCTION subprogranm. 
None of the dummy arguments may appear in an EQUIVALENCE or COMMON 
statement. 


Type Specification of the FUNCTION Subprogram 


In addition to implicitly declaring the type of a FUNCTION name by 
the predefined convention, there exists the option of explicitly 
specifying the type of a FUNCTION name within the FUNCTION statement. 
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However, the type of a FUNCTION name may not be declared in an Explicit 
Specification statement in that subprogram. 


NE ANY AE ES TE EE RS EO OT AS EN COS aN AS ED ERAS GREP ILE GE CEE GREE ED RED REE UR AY SREY ARSED RID FEY ANY GRP OD TED GREY EEE GE <ED AONED ERND GE SD CIES OE CUEY AO emer ceser eee CAE tees CUR EP SRS eee eee ee CTE ees Ee ee a oe GE ED ames ED SEE ae eR 


P 
| General Form 


Type FUNCTION name (a1, Aa,--- san) 
Where: Type is INTEGER, REAL, Or DOUBLE PRECISION. 
name is the name of the FUNCTION subprogram. 
AgeGaeeee,an are nonsubscripted variable or array names, or 


the dummy names of a SUBROUTINE or another FUNCTION subpro- 
gram. (There must be at least one argument.) 


[eo ae ee ee. ae en oe ee ee ae 


Example 1: 


REAL FUNCTION SOMEF (A,B) 


SOMEF = A*¥*2 + B*¥*2 


RETURN 
END 


Example 2: 
INTEGER FUNCTION CALC(X,Y,2Z) 


CALC = Xt+Y¥+Z**2 


«ee 


RETURN 
END 


Explanation: 


The FUNCTION subprograms SOMEF and CALC in Examples 1 and 2 are 
declared as type REAL and INTEGER, respectively. 


RETURN and END Statements in a Function Subprogram 


All FUNCTION subprograms must contain an END statement and at least 
one RETURN statement. The END statement specifies, for the compiler, 
the physical end of the subprogram; the RETURN statement signifies a 
logical conclusion of the computation and returns any computed value and 
control to the calling program. There may, in fact, be more than one 
RETURN statement in a FORTRAN subprogram. 


72. 


Example: 


FUNCTION DAV (D,E,F) 
IF (D-E) 10, 20, 30 
10 A = Dt2.0*E 


5 A = F+t2.0*E 


20 DAV = At+Bt*2 


RETURN 
30 DAV = B¥¥2 


RETURN 
END 


SUBROUTINE SUBPROGRAMS 


The SUBROUTINE subprogram is Similar to the FUNCTION subprogram in 
many respects: the rules for naming FUNCTION and SUBROUTINE subprograms 
are the same, they both require an END statement, and they both contain 
the same sort of dummy arguments. Like the FUNCTION subprogram, the 
SUBROUTINE subprogram is ae_eset of commonly used computations, but it 
need not return any results to the calling program, as does the FUNCTION 
subprogram. 


The SUBROUTINE subprogram is called by the CALL statement, which 
consists of the word CALL followed by the name of the subprogram and its 
parenthesized arguments. 


| General Form 


SUBROUTINE name (a4¢a2 «Ase oe °*an) 


RETURN 
END 


where: name is the subprogram name (see the section "Naming 
Subprograms"). 





A4742ne¢43eeeeyan are nonsubscripted variable or array names, 
or the dummy names of other SUBROUTINE or FUNCTION subpro- 
grams. (There need not be any arguments.) 


Ee Ee ED EE OREIEY ETE AE ED SISED AEEEED <EEED ER CEE mE GEE EER SERED qu QTE eee CENTRE EEE CHOED coe GIES thE men exgEe GEE) GERD CREED <eREE SEED GEE GREP cee <TD CED GD GeEUD GED cueee denne ee ED Ce ee cee Oe ee eee oe ERD ee ED DEED ED eee come “SED que ques Ee <tmee quine ce aes ee ae 


ee come cere eee sees Se eee eee TE ee ee ee ee Ee ee a 


Since the SUBROUTINE is a separate subprogram, the variables and 
statement numbers within it do not relate to any other program. 


The SUBROUTINE subprogram may contain any FORTRAN statement except a 
FUNCTION statement or another SUBROUTINE statement. 
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The SUBROUTINE subprogram may use one or more of its arguments’ to 
return values to the calling program. Any arguments so used must appear 
on the left side of an arithmetic statement or in an input list within 
the subprogram, as arguments of a CALL statement or as arguments in a 
function reference. The SUBROUTINE name must not appear in any other 
Statement in the SUBROUTINE subprogranm. 


The arguments (a4, @as A3,-+-2,an) May be considered dummy variable 
names that are replaced at the time of execution by the actual arguments 
supplied in the CALL statement. The actual arguments must correspond in 
number, order, and type to the dummy arguments. The array size must 
also be the same. Dummy arguments may not appear in an EQUIVALENCE or 
COMMON statement within the subprogram. 


Example: The relationship between variable names used as arguments in 
the calling program and the dummy variable used as arguments in the 
SUBROUTINE subprogram is illustrated in the following example. The 
object of the subprogram is to “copy" one array directly into another. 


Calling Program SUBROUTINE Subprogram 
DIMENSION X(100),Y(100) 
‘ SUBROUTINE COPY(A,B,N) 
‘e DIMENSION A (100) ,B(100) 
; DpDO 10 T= 1, N 
CALL COPY (X,Y,K) 10 BCI) = A (TI) 
‘ RETURN 
‘ END 


@ 


CALL Statement 


The CALL statement is used to call a SUBROUTINE subprogram. 


}---------------------------------------------------------------------- i 
| | 
| CALL name (a4,82,a83,--+-+,aAn) | 
| | 
| Where: name is the name of a subroutine subprogram. | 
| | 
| Q4,42sa3geeee,7an are the actual arguments that are being | 
| supplied to the subroutine subprogram. | 
i eh uch a 8 a a ee ate J 
Examples: 


CALL OUT 

CALL MATMPY (X,5,40,Y, 7,2) 
CALL QDRTIC (X,Y,Z,ROOT1, ROOT2) 
CALL SUB1(X+Y¥*5, ABDF, SINE) 
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The CALL statement transfers control to the subroutine subprogram and 
replaces the dummy variables with the value of the actual arguments that 
appear in the CALL statement. The arguments in a CALL statement may be 
any Of the following: 


1. any type of constant 

- any type of subscripted or nonsubscripted variable 
- an array name 

- an arithmetic expression 

- the name of a FUNCTION or SUBROUTINE subprogram 


OI & WwW Nh 


The arguments in a CALL statement must agree in number, order, and 
type with the corresponding arguments in the subroutine subprogram. The 
array sizes must also be the same in the subroutine and the calling 
programs. If an actual argument corresponds to a dummy argument that is 
defined or redefined in the subprogram, the actual argument must be a 
variable name, subscripted variable name, or array name. 


A subprogram cannot define dummy arguments when the subprogram 
reference causes those arguments to be associated with other dummy 
arguments within the subprogram or with variables in COMMON. For 
example, if the external function DERIV is defined as 


FUNCTION DERIV (X,Y, 2) 
COMMON W 


and if the following statements are included in the calling source 
program 


COMMON B 


C = DERIV (A,B, A) 


then xX, Y, Z, and W cannot be defined (1i.e., cannot appear to the left 
of an equal sign in an arithmetic statement) in the function DERIV 
because the actual argument list causes both A and B to be associated 
with more than one value. 


RETURN Statement in a SUBROUTINE Subprogram 


Gs STS aa a aa aaa 1 
| General Form | 
fa a nn nr er a { 
| 
| RETURN | 
a a a a a ) 


This is the exit from a subprogram. The RETURN statement signifies 
the conclusion of a computation and returns control, and any values 
requested, to the calling program. In aA main program, a RETURN 
Statement performs the same function as a STOP statement. There may be 
several RETURN statements in a subprogram. 
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EXTERNAL Statement 


A aa a a a aaa ac a a | 
| General Form | 
fem a { 
| | 
| EXTERNAL a,b,c,... | 
| : | 
| Where: a,b,c,.-. are names of subprograms that are passed as | 
| arguments to other subprograms. | 
Masses i a at saa a ee J 


The EXTERNAL statement is a specification statement and must appear 
prior to any executable statement in the source program. 


The name of any subprogram that is passed as an argument to another 
subprogram must appear in an EXTERNAL statement in the calling program. 
For example, assume that SUB and MULT are subprogram names in the 
following statements: 


Example 1: 
Calling Program Subprogram 
ie SUBROUTINE SUB(K,Y,2Z) 
P IF (K) 4,6,6 


° 4 D= Y¥ (K,2Z**2) 
EXTERNAL MULT 


CALL SUB (J, MULT,C) 6 RETURN 
° END 


@ 


Explanation: 


In this example, the subprogram name MULT is used as an argument in 
the subprogram SUB. The subprogram name MULT is passed to the dummy 
variable Yas are the variables J and C passed to the dummy variables K 
and Z, respectively. The subprogram MULT is called and executed only if 
the value of K is negative. 


Example 2: 


CALL SUB (A,B,MULT (C,D),37) 


e@ 
@ 


@ 


Explanation: 


In this example, an EXTERNAL statement is not required because the 
subprogram named MULT is not an argument; it is executed first and the 
result becomes the argument. 
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APPENDIX A: SOURCE PROGRAM CHARACTERS 


Source programs may be coded in either BCDIC or EBCDIC’ character 
codes. Mixing of the two, however, is not allowed. 


EBCDIC EBCDIC 

or BCDIC or BCDIC 
Alphabetic Card Numeric Card 
Characters Punches Characters Punches 








42-1 
L2e2 
12-3 
12-4 
12-5 
12-6 
Lard 
12-8 


WAN AMNLWNHF O 
ODMDNINMNLWHF SO 


UVNK MSE GAHAN DWONOBABEAHRYUH MAA tOAOW Pp 
my 
! be 
WO DOYIA UO FW ! 
wn 


iy ora ck des Ge da 


1-8-3 


Special EBCDIC Card BCDIC Card Dual 
Characters Punches Punches Characters 





(blank) (no punch) (no punch) 

- 11 11 

> 11-8-4 11-8-4 
12-8-3 12-8-3 

(comma ) 0-8-3 0-8-3 
12-8-6 12 & (ampersand) 
0-1 0-1 
8-6 8-3 
12-8-5 0-8-4 
11-8-5 12-8-4 

(apostrophe) 8-5 8-4 


wANN +s 


| 


& TT 3a + 
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APPENDIX B: BASIC FORTRAN IV IMPLEMENTATION DIFFERENCES 


The differences among the four implementations of the Basic FORTRAN 
IV language are minor except for the absence of the direct access 
input/output statements in Basic Programming Support FORTRAN IV. These 
differences are indicated in the body of this publication by footnotes 
and, in addition, are summarized in Table 5. 


The following abbreviations are used in Table 5: 


OS - Operating System FORTRAN IV (E) 

DOS - Disk Operating System FORTRAN IV 

TOS Tape Operating System FORTRAN IV 

BPS Basic Programming Support Tape System FORTRAN IV 


78 


Table 5. 


Pen ence ete ere Sem 
| 


Direct-access 
[f input/output 
istatements 


|Key word 
fand blank 
|restriction 


{Maximum array 
{size 


+ a ec le een ee rm a a 


| Maximum 
|record size 


{Order of 
| Specification 
|statements 


Tape ene OEE se: eee eet eee eee OEE gee Se eee oe 


| Direct-—access 


TM Ge en ee eee —~T 


Implementation Differences 
peered eee 


| DOS/TOS [BPS 


|Direct-access I/0|Direct-access I/0O| 


JI/O is available.|is available. A | 
| {program with | 
| |direct-access I/O| 
i Jmay be compiled | 
| [using either DOS | 
| jor TOS, but it | 
| |must be executed | 
| Jusing DOS. i 
}~--~----------~--}----------------- + 
{Control card op- |The restriction 
Jtion may be used |has been removed; | 
[to either keep or[there is no op- | 
[remove the re- Jtion to retain | 


|striction. See |{it. | 
jthe programmer's | | 
| guide. [ | 
-}--------------—-- }-~---~-----------+ 
}131,068 storage |32,767 storage | 
| locations | locations | 
}~---------------- }~----~----------- + 


{Depends upon {255 characters | 
[the input/output |per record | 
| device in use; | 

|See the program- 

{mer's guide. = 


| Specification 
|statements that 
{describe data 
[must precede any 
|statements which 
|xrefer to that 


| Specification | 
{statements that | 
{describe data | 
[must precede any | 
| statements which | 
{refer to that | 


|data. EQUIVALENCE|data. EQUIVALENCE|ments. 


[statements must [statements must | 
{follow the ex- | follow the ex- | 
[plicit specifica-|plicit specifica~| 
[tion, COMMON, and|tion, COMMON, and| 
|DIMENSION state- |DIMENSION state- | 
[ments but do not |ments but do not | 
jhave to follow [have to follow [ 
|lthe DEFINE FILE |the DEFINE FILE | 
j|statement. | statement. [ 


EP epee come ee eaee same qme> ent anee SEnDe cEEe qnED ce ED eae een cD er ERE cee <eee EE ED comes QED oe GY ED qe cee eee ES ND EEE KEEP EDD ES ED ENED ENDS ED FEED RS ED ED EN AED ee I I a ed 


is not available. | 
| 
| 
| 
| 
| 
| 


|The restriction 


| 
has been removed; | 
there is no op- | 
tion to retain | 
| 
| 
| 


32,767 storage { 
locations 
Ss Bs ep ed mee eran eye J 
255 characters 
per record 


ee Re ee eee anne Soe aD RD eee eee Ue RD eeD oEe  er 


The EQUIVALENCE 
Statements must 
follow the ex- 
plicit specifica- 
tion, COMMON, and 
DIMENSION state- 


sche nin mimi lori mip STO gnc cits” SNe ces “mm is mis mtn wang, calcite; aioe Gemmm outline 


{Subprogram names|Any valid FORTRAN|Any valid FORTRAN|Any valid FORTRAN| 


{2 sae See ee eee ee 


Appendix B: 


js as asa ak ls ee eel er eer oes 


|mame may be used |name may be used. | 
[unless key word 
{restriction is 
jretained; key 
{words may not 
jthen be used. 


a ee ee ee 


SUP ee RD ED AES ED Ee cH ERNE ee ee TD <e em ae 


Basic FORTRAN IV Implementatio 


name may be used. 


| 
| 
| 
| 
| 
| 
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APPENDIX C: FORTRAN SUPPLIED SUBPROGRAMS 


The FORTRAN supplied subprograms are either in-line or out-of-line. 
An in-line subprogram is inserted by the FORTRAN compiler at any point 
in the program where the subprogram is referenced. The in-line 
subprograms are mathematical function subprograms. These subprograms 
are listed in Table 6. 


The out-of-line subprograms are located on a library. These subpro- 
grams are mathematical function subprograms and service subprograms. 
The out-of-line mathematical function subprograms are listed in Table 7; 
out-of-line service subprograms are listed in Table 8. A detailed 
description of all out-of-line subprograms is contained in the publica- 


tion IBM System/360 Operating System: FORTRAN IV (E), Library Subpro- 
grams. 


Note: Variables used as arguments of any mathematical function subpro- 
gram must be defined in accordance with the function in which they 
appear. This definition is accomplished either with the Explicit 
specification statement or with the predefined convention. 





Table 6. In-Line Mathematical Function Subprograms 


iii aaa an Canna amin GEER aca al aioe ae a ea eae Rane: (Gai ne Tae is: 1S ea aaa a 
JEntry | |No. of| Argument |Function | 
| Function [Name |Definition {| Arg. | Type jValue Type| 
tt ee ee 
{Absolute value J IABS ||Arg| | 1 | Integer | Integer | 
i | ABS | | 1 |Real | Real | 
{ {DABS | } 1 | Double | Double | 
| | | | | Precision| Precision| 
~--~-------------}------}----------------} ------ f----------}---------- 
| Float | FLOAT |Convert from {| 1 j Integer | Real [ 
| {|DFLOAT|integer to real | 1 | Integer | Double [ 
i i | | { {Precision | 
~---------~------}------}----—-----—----} —---- }----------}----------{ 
| Fix {IFIX |Convert from j 1 | Real | Integer { 
| | {real to integer | | | | 
~------~---------}------}---------------- + ---- f---------- f= -- 
[Transfer of sign |SIGN |Sign of Arga |} 2 |Real |Real | 
| | [times |Arg, | | | | 
; | ISIGN | | 2 | Integer | Integer [ 
i |DSIGN | | 2 | Double | Double | 
| | | | | Precision| Precision| 
}-----------------}------}----------------}------ f-------------------- 
| Positive | DIM JArg,-Min(Arg,, | 2 | Real | Real 
| difference [| IDIM | Arga) | 2 | Integer | Integer [ 
mannan =f fn nnn fn nnn nnn nnn nnn nnn 
{Obtaining most |SNGL | | 1 [Double [Real | 
|significant part | { [ | Precision| | 
jof a Double | [ i | | 
| Precision | | i i | [ 
| argument | l ff | 
}-----------------}------}----------------4 —---- f----------}----------] 
|Express a Real [DBLE | | 1 | Real | Double | 
jargument in | [ | Precision| 
[Double Precision | l [ | | | 
| form | | | | | | 
Eee eae a ee ene B Epeeenee-arestonrne tee eS p ietaes ee eae dE ee een ey Lennar mee! 


Table 7. Out-of-Line Mathematical Function Subprograms 


aaa earn a a, ae mire Mtge ey nee eC et ee ee er eee 1 
| [Entry | | No. of | Argument [Function 
| Function [Name |Definition | cae i Type |Value Type| 
ep ee ee eee oN See ae sar Sa ea 1 se Ra cine SRR eames" Gaui aimee cia 
| Exponential | EXP jeArd ; 4 {Real |Real | 
| J|DEXP |eatg | 1 | Double | Double | 
| | | | Precision| Precision| 
~---~--~----~----}----~-}---—------------ }----—— f---------- f----------4 
|Natural Logarithm|ALOG | Ln(Arg) { 1 |Real {Real | 
{ {|DLOG |{Ln(Arg) ij 1 {Double | Double | 
| | | | | Precision| Precision| 
}-----------~-----}~-----}-----~-------~--~------}----------}----------4 
[Common Logarithm |ALOG10|40g,, (Arg) j} 1 | Real | Real | 
i | DLOG10|Lo (Arg) ; 1 |Double |Double { 
Gao g e ° s . 
{ | | i { Precision| Precision| 
~----------------}------ }--~--~-~--~---~-}------ }----------}----------+ 
{| Arctangent |ATAN [arctan (Arg) j| 1 | Real | Real | 
| |DATAN [arctan (Arg) {| 1 {Double | Double | 
{ | | | | Precision| Precision| 
}---~----------—--}------}---------------- }------ +---------- +---------- { 
| Trigonometric | SIN [sin (Arg) j 1 |Real |Real i 
j|Sine |DSIN |sin(Arg) } 1 | Double | Double | 
| | | Precision| Precision| 
}--------------—--}—---- }-~----------~--- }----~-}----------}---------- { 
| Trigonometric | COS | cos (Arg) j 1 | Real | Real | 
| Cosine {DCOS [cos (Arg) { 1 {Double | Double | 
| | | i | Precision| Precision| 
~~ -----~ == == ===} =f nnn nf nnn { 
{Square Root {SORT | (Arg) ; 1 | Real | Real | 
r |DSORT | (Arg)”? 1 |Double | Double l 


| Precision] Precision] 


t t + t 
| Hyperbolic |TANH [tanh (Arg) 1 |Real {Real 
| Tangent J|DTANH {tanh (Arg) 1 | Double | Double | 
| | { | Precision| Precision| 
~--~---~-~--~-~~-}~----— f-----=-------- =} -----~ f---------- f---------- 
|Modular | MOD [Arg, (mod Arga) | 2 | Integer | Integer | 
j Arithmetic JAMOD | { 2 |Real | Real | 
| (Remaindering) [|DMOD | {| 2 {Double | Double | 
| | | i | Precision| Precision| 
Sree Se wae oe ed ee See ee ee a ee ee 
| Truncation { INT [Sign of Arg | 1 {Real [Integer | 
| { Jjtimes largest [ | [ i 
| [integer <|Arg| | | | 
| [|AINT | { 1 | Real {Real | 
| {| IDINT | {| 1 | Double | Integer | 
| ] | i | Precision| | 
-~---------------- }------}-~-~--------—----}------}----------}---------- 
| Largest value [AMAXO {Max (Arg,, | 22 |Integer [Real | 
{| AMAX1 | aeezsAKGyn) | 22 #4«|Real {Real | 
| |MAXO | | 22 =|iInteger | Integer | 
H |MAX1 | {| 22 |Real | Integer | 
| | DMAX1 | j 22 #.|Double | Double | 
i | | ( | Precision| Precision| 
~~~--------------}------ }-----~----------}~-----}----------}---------- 
|Smallest value [AMINO |Min (Arg,, | 22 |Integer {Real | 
| | AMIN1 | eo2,AYgn) | 22 [Real | Real | 
| {MINO | | 22 |Integer | Integer | 
| [MINI | { 22 {Real | Integer i 
| {DMIN1 | | 22 [Rouble | Double [ 
| | | i | PEee oo: poreee ron 
a ie Baretta a rea ee ee se 
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Table 8. Out-of-Line Service Subprograms 


oe ee ea ee ee ba aR we oe, Be 


| Function }Call Statement 


AS END ED EN TED OED ED ED ED EE ED ee LD SRD EEE RR RENE UD ETS ED I ED CE ED EE ee Ee ee ee ee 


[Alter status of|CALL SLITE(i) 
[sense lights | 


cfeee exces SRD SUN UO ee ae 


{Test and record|CALL SLITET (i, j) 


{status of sense] 
{lights | 


[Dump storage on|CALL DUMP (a,,b,,f., 


jthe output data| 
jset and termin-] 
Jate execution 


j|Dump storage on|CALL PDUMP (a,,by,f,, 


jthe output data] 
Jset and contin-| 
Jue execution | 


°°*7anrDn, fn) 


2+228aneDnsfn) 


sn ge enh tt i's im i Sin si si ini ner inns She ni si ih si Sen 1 


+ 
{Argument Information | 


| 

| 

| 

|If i+ 0, the four sense| 
{| lights are turned off. | 
{If i= 1,2,3,4, the corres-| 
{| ponding sense light is | 
{| turned on. | 
ee eee 
{i is an integer expression | 
| that has a value of 1, 2, 3,| 
| or 4 and indicates which| 
| sense light to test. 

{3 is an integer variable that| 
{| ais set to 1 if the sense 
| light was on; or to 2 if the| 
| sense light was off. | 
4 


a and b are variables that 
indicate the limits of stor-| 
age to be dumped. (Either a| 
or b may be the upper or| 
lower limits of storage, but| 
both must be in the same| 
program or subprogram or in| 
COMMON. ) 

f indicates the dump format| 
and may be one of the fol- 


lowing: 
4 - integer 
5 - real 
6 - double precision 
a, b, and £ are as previously 
defined f 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| | 
| | | 
| 0 - hexadecimal | 
| | 
| | 
| ( 
4 
| 
| 
| 
| 


| 
or DUMP. | 
| 
| 


sc aa tesa i em ae a -——-—-—-——-—-— ———- -—- —-—-—- —- -—- -—- —- + -- - +--+ + 


[Test for divide|CALL DVCHK (4) 
|check exception| 


|Test for expon-|CALL OVERFL(j) 
[ent overflow or| 
junderflow | 


— 
i 
i 
| 
i 
! 
| 
| 
i 
i 
{ 
| 
{ 
| 
! 
! 
i 
! 
| 
{ 
! 
! 
| 
i 
| 
| 
| 
| 
| 
i 
| 
| 
| 
| 
! 


| Terminate exec-]CALL EXIT 
J ution | 


if 
' 
| 
1 
| 
! 
| 
| 
i 
| 
i 
! 
| 
l 
J 
- 
| 
i 
! 
i 
| 
! 
1 
i 
| 
{ 
| 
! 
| 
| 
| 
| 
I 
! 
{ 
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ek 
| 
i 
| 
| 
| 
! 
| 
| 
{ 
| 
i 
| 
! 
| 
| 
{ 
i 
| 
i 
| 
| 
| 
| 
| 
| 
| 
i 
| 
| 
amis 


{7 is an integer variable that | 
is set to 1 if the divide-| 
check indicator was on; or| 
to 2 if the indicator was| 
off. After testing, the | 
divide-check indicator is | 
turned off. | 

j is an integer variable that | 
is set to 1 if an overflow | 
condition exists; to 2 if no| 
overflow condition exists; | 
or to 3 if an underflow 
condition exists. After | 
testing, the overflow indi-| 
cator is turned off. | 
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APPENDIX D: SAMPLE PROGRAMS 


SAMPLE PROGRAM 1 


The sample program (Figure 2) is designed to find all of the prime 
numbers between 1 and 1000. A prime number is an integer that cannot be 
evenly divided by any integer except itself and 1. Thus 1, 2, 3, 5, 7, 
11,... are prime numbers. The number 9, for example, is not a prime 
number since it can evenly be divided by 3. 
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Figure 2. Sample Program 1 
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SAMPLE PRCGRAM 2 
The n points (x;, yj) are to be used to fit an m degree polynomial by 
the least-squares method. 


Y = Agt a4qxX + AQX* + 22oe + Ayx™M 


In order to obtain the coefficients ag, a4,«+-, am, it is necessary to 
solve the normal equations: 


(2) W1ag + Waaq + 22. + Wyeaam = Za 
(m+ 1) Wma + Wm 4124 + oe @ + Wam@m = 4m 
where: 
n 
No > a a ee Ve 
i=1 
n n 
a ee Zy =X YiXj 
n n 
Wa = 2% xi" Zo =X Y;xj? 
x n 
; 2m = & yyxi™ 
° i=1 
n 
Wam = > x., 40 
i=1 


After the W's and Z's have been computed, the normal equations are 
solved by the method of elimination which is illustrated by the 
following solution of the normal equations for a second degree polynomi- 
al (m = 2). 


(1) Wo@o + Waad + Woag = Z, 
(2) W4aoq + Waa, + Waae = Zy 
(3) W2a, + Wsa, + Wia2 = Z2 


The forward solution is as follows: 
1. Divide equation (1) by Wy 


2- Multiply the equation resulting from step 1 by W, and subtract from 
equation (2). 


3. Multiply the equation resulting from step 1 by Ws. and subtract from 
equation (3). 
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The resuiting equations are: 


(4) @g + banaq + Dynae = Day 
(5) B22aq + baga2 = Day 
(6) Bs2aq + ba3a2 = Day 
where: 
bea = WA, ba3 = Wo/Wo, Dan = Zo/Wo 


i 


Da2 = Wa-b12Wa4 , bas W3-bDaaWa , Day = Za-Da, Wy 
D32 = Wa-ba2Wo , Daz = Wa-bDaaWe , Day = Za-biy We 


Steps 1 and 2 are repeated using equations (5) and (6), with bj, and bao 
instead of Wy. and W,. The resulting equations are: 


(7) aq + C2382 = C2y 
(8) C33a€2 = C34 
where: 
C23 = bas/bDaa s Cag = Day/bDaa 
C33 = ba3-Ca3D32 , Cau = D3y-CaqgbD32 


The backward solution is as follows: 


(9) Q@2 = C3a4/C33 | from equation (8) 
(10) Aq = Cay~C23a2 from equation (7) 
(11) a = Day -b428q-Da 382 from equation (4) 


Figure 3 is a possible FORTRAN program for carrying out the 
calculations for the case: n = 100, m < 10. Wo, Wa, Wary cee, Wam are 
stored in W(1), W(2), W(3), «-., W(2M+1), respectively. Zoe Lae Gaze 
eco, Zp are Stored in Z (1), Z (2), 2 (3), «--, Z(Mt1), respectively. 
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Figure 3. Sample Program 2 (Continued) 


The elements of the W array, except W(1), are set equal to zero. 
W(1) is set equal to N. For each value of I, Xj and Yj are _ selected. 
The powers of Xj are computed and accumulated in the correct W counters. 
The powers of Xq4 are multiplied by Y; and the products are accumulated 
in the correct Z counters. In order to save machine time when the 
object program is being run, the previously computed power of Xj is used 
when computing the next power of Xj. Note the use of variables as index 
parameters. By the time control has passed to statement 17, the 
counters have been set as follows: 


N 
W(1) = N Z(1) =Z% Yy 
I=1 
N N 
W(2) =X Xr Z2(2) =X YyzXy 
I=1 I=1 
N N 
w(3) =X Xr? Z(3) = XZ YyXyz? 
I=1 I=1 
a N o 
‘ Z(M+1) = YyX, 
- I=1 
N 
W(2M#1) =X x; ™ 
I=1 
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By the time control has passed to statement 23, the values Of Wo, W,, 
eces Wamts1 have been placed in the storage locations corresponding to 
columns 1 through M + 1, rows 1 through M + 1, of the B array, and the 
values Of Z , Za, eee, Zy have been stored in the locations correspond- 
ing to the column of the B array. For example, for the illustrative 
problem (M = 2), columns 1 through 4, rows 1 through 3, of the B_ array 
would be set to the following computed values: 


Wo W, Wo ae 
Wy We W3 Ls 
Wo W3 Wa Za 


This matrix represents equations (1), (2), and (3), the normal 
equations for M = 2. 


The forward solution, which results in equations (4), (7), and (8) in 
the illustrative problem, is carried out by statements 23 through 31. 
By the time control has passed to statement 33, the coefficients of the 
AI terms in the M + 1 equations which would be obtained in hand 
calculations have replaced the contents of the locations corresponding 
to columns 1 through M+1, rows 1 through Mt1, of the B array, and the 
constants on the right-hand side of the equations have replaced the 
contents of the locations corresponding to column Mt+2, rows 1 through 
M+t1, of the B array. For the illustrative problem, columns 1 through 4, 
rows 1 through 3, of the B array would be set to the following computed 
values: : 


1 baz bag Dow 
0 1 C23 Coa 
0 0 C33 Cau 


- This matrix represents equations (4), (7), and (8). 


The backward solution, which results in equations (9), (10), and (11) 
in the illustrative problem, is carried out by statements 33 through 40. 
By the time control has passed to statement 41, which prints the values 
of the A9 terms, the values of the (M+1)#A; terms have been stored in 
the M + 1 locations for the A array. For the illustrative problem, the 
A array would contain the following computed values for ag, a4, and @g, 
respectively: 


Location Contents 

A (3) C3au/C33 

A (2) C24-C2342 

A (1) Dea ~Dq 2814-Dq 382 


The resulting values of the Ay terms are then printed according to 
the FORMAT specification in statement 2. 
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APPENDIX E: FORTRAN IV FEATURES NOT IN BASIC FORTRAN IV 


The following statements and features in FORTRAN IV are not in Basic 
FORTRAN IV: 


ASSIGN 

BLOCK DATA 

Labeled COMMON 

COMPLEX 

DATA 

More than three dimensions 
Adjustable dimensions 

Assigned GO TO 

Logical IF 

LOGICAL 

PRINT b,list 

PUNCH b,list 

READ b, list 

END and ERR parameters in a READ 
Generalized Type statement (But note that DOUBLE PRECISION is 
provided as an explicit type) 
IMPLICIT 

Call by name 

Literal as argument of CALL 

ENTRY 

RETURN i (i not a blank) 
NAMELIST 

PAUSE with literal 

G, L, and Z format codes 
Complex, logical, literal, and hexadecimal constants 
Generalized subscript form 
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A format code 44-46 
ABS 78 
Addition 
(see Arithmetic operators) 
AINT 81 


ALOG 81 

ALOG10 81 
AMAXO 81 
AMAX1 81 
AMINO 81 
AMIN1 81 
AMOD 81 


Arithmetic expressions 14-18 
mode of 14,15-16 
rules for writing 15-18 
Arithmetic IF 24-25 
Arithmetic operators 14-16 
Arithmetic statement 7,22 
Arrangement of arrays in storage 21 
Arrays 18-21 
arrangement in storage 21 
Size declaration 21 
Subscripted variables 19 
Subscripts 20-21 
ASA basic FORTRAN 5 
ATAN 81 


BACKSPACE statement 52 
BaSic FORTRAN IV 5-6,78-79 
Blank fields 

(see X format code) 
Blank lines 

(see Carriage control) 
Blanks 8-9 


C (see Comments line) 
CALL statement 74-75 
Carriage control 51 
Character card punch codes 77 
Coding form 8 
Comments line 8 
COMMON statement 62-64 
Compilers 5 
Computed GO TO statement 24 
Constants 9-11 
double precision i11 
integer 9 
real 10 
Continuation lines 7 
CONTINUE statement 28-29 
Control statements 7,23-30 
arithmetic IF 24-25 
computed GO TO 24 
CONTINUE 28-29 
DO 25-28 
END 30 
PAUSE 29 
STOP 30 
Conversion codes (see Format codes) 
Cos 81 | 


INDEX 


D format code 41,43 
DABS’ 80 
Data set (defined) 31 
Data set reference number (defined) 31 
DATAN 81 
DBLE 80 
DCOS 81 
Declaring the size of an array 21 
DEFINE FILE statement 53-55 
DEXP 81 
DFLOAT 80 
Differences among implementations 78-79 
DIM 80 
DIMENSION statement 61 
Division 
(see Arithmetic operators) 
Direct-access I/O statements 53-59 
DEFINE FILE 53-55 
FIND 58 
READ 55-56 
WRITE 56-58 
DLOG 81 
DLOG10 81 
DMAX1 81 
DMINi 81 
DMOD 81 
DO loops 25-28 
DO statement 25-28 
increment 25,26,27,28 
initial value 25,26,27,28 
programming considerations 27-28 
range 25,27-28 
test value 25,26,27,28 
variable 25,26,27,28 
Double precision constants 11 
DOUBLE PRECISION statement 60-61 
DSIGN 80 
DSIN 81 
DSORT 81 
DTANH 81 
DUMP 82 
DVCHK 82 


E format code 41,43 
END FILE statement 52 
END statement 30 

in FUNCTION subprograms 72-73 
EQUIVALENCE statement 64-66 
EXIT 82 
EXP 81 
Explicit specification statement 60-61 
Exponentiation 

(see Arithmetic operators) 
Exponents 

(see E and D format codes) 
Expressions 14-18 
EXTERNAL statement 76 


F format code 41,42-43 


Index 91 


FIND statement 58 


FLOAT 80 
Format codes 41-51 
A code 44 


carriage control 51 
D and E codes 41,43 
F code 41,42-43 
H code 47 
I code 41-42 
numeric codes 41-43 
scale factor - P 49-51 
T code 48-49 
X code 48 
FORMAT statement 37-51 
format codes 41-51 
FORTRAN record 38-40 
literal data 46-47 
use of 38-41 
FORTRAN supplied subprograms 67,80-82 
EXIT, DUMP, and PDUMP 82 
in-line mathematical function 
Subprograms 81 
machine indicator tests 82 
out-of-line mathematical function 
Subprograms 82 
FUNCTION subprograms 67,/70-73,80-81 
FORTRAN supplied 80-82 
type specification of 67,71-72 
use of RETURN and END 72-73 
Functions 68-73 
definition 68 
FUNCTION subprograms 70-73 
reference to 68 
Statement functions 68-69 
type specification of 67,71-72 


GO TO statements 23-24 
computed 24 
unconditional 23 


H format code 47 
Hierarchy of operations in an arithmetic 
Statement 16-17 


I format code 41-42 


IABS 80 
IDIM 80 
IDINT 81 
IFIX 80 


IF statement 24-25 
Implementation differences 78-79 
Input/output statements 7, 31-59 

BACKSPACE 52 

DEFINE FILE 53-55 

direct access 53-59 

END FILE 52 


FIND 58 
READ 32-35,55-56 
REWIND 52 


sequential 31-37,52 
WRITE 35-37,56-58 
INT 81 
Integer constants 9 
Integer division 18 
INTEGER statement 60-61 
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I/O lists 32-37,54-55,57 
ISIGN 80 


Library 
(see FORTRAN supplied subprograms) 
Literal data 46-47 


Machine indicator tests 82 
Mathematical function subprograms 80-81 
MAXO 81 


MAX1 81 

MINO 81 

MINi 81 
Mixed mode 6 
MOD 78,81 


Mode of an arithmetic expression 14,15-16 
Multiplication 
(see Arithmetic operators) 


Naming subprograms 67 
Nesting DO loops 27-28 
Numeric characters 77 
Numeric format codes 41-43 

D format coae 43 

E format code 43 

F format code 42-43 

I format code 41-42 


Operands 15 
Operators 
(see Arithmetic operators) 
Order of computation in an arithmetic 
expression 16-17 
OVERFL 82 


P format code 49-51 

PAUSE statement 27 

PDUMP 82 

Predefined specification 
(convention) 13 

Programming considerations 
in DO loops 27-28 
using direct access 54,58-59 


Range of a DO statement 25,27-28 
READ statement 32-35,55-56 
READ (a) list 33-34 
READ (a,b) list 32-33 
READ (a'r,b) list 55-56 
Real constant 10 
REAL statement 60-61 
RETURN statement 72-73,75 
REWIND statement 52 


Sample programs 83-87 
Scale factor - P 49-50 
Sense light subprograms 82 
Sequential I/O statements 31-37,52 
BACKSPACE 52 
END FILE 52 
READ 32-35 


REWIND 52 
WRITE 35-37 
SIGN 80 
SIN 81 


SLITE 82 
SLITET 82 
SNGLE 80 
Source program characters 77 
Special characters 77 
Specification statements 7,60-66 
COMMON 62-64 
DIMENSION 61 
EQUIVALENCE 64-66 
Explicit specification 60-61 
SQRT 81 
Statement functions 67,68-69 
STOP statement 30 
Subprograms 67-76, 80-82 
FORTRAN supplied 80-82 
FUNCTION subprogram 70-73,80-81 
functions 68-73 
naming 67 
Statement functions 68-69 
SUBROUTINE subprograms 73-76,82 
SUBROUTINE subprograms 67, 73-76, 82 
Subscript 18,20-21 


Subscripted variables 18,20-21 
Subtraction 
(see Arithmetic operators) 


T format code 48-49 

TANH 81 

Type declaration 13 

Type specification 60-61,71-72 


Unconditional GO TO 23 


Variables 12-13,18-21 
Subscripted 18-21 
type specification 13 
variable names 12 


WRITE statement 35-37,56-58 
WRITE (a) list 37 
WRITE (a,b) list 36 
WRITE (a‘'r,b) list 56-58 


X format code 48 





Index 
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